Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 27 October 201 1 (MN KT^ style file v2.2) 



Can galactic outflows explain the properties of Lya emitters? 



Alvaro Orsi^'^'^*, Cedric G. Lacey^ and Carlton M. Baugh^ 

1. Departamento de Astronomta y Astrofisica, Pontificia Universidad Catolica, Av. Vicuna Mackenna 4860, Santiago, Chile. 

2. Centra de Astro-Ingenieria, Pontificia Universidad Catolica, Av. Vicuna Mackenna 4860, Santiago, Chile. 

3. Institute for Computational Cosmology, Department of Physics, University of Durham, South Road, Durham DHl 3LE, UK. 



o 

(N 
» 

o 
O 

(N 

o 

u 

6 



o 
d 



27 October 2011 



ABSTRACT 

We study the properties of Lya emitters in a cosmological framework by computing the es- 
cape of Lya photons through galactic outflows. We combine the GALFORM semi-analytical 
model of galaxy formation with a Monte Carlo Lya radiative transfer code. The properties 
of Lya emitters at < z < 7 are predicted using two outflow geometries: a Shell of neutral 
gas and a Wind ejecting material, both expanding at constant velocity. We characterise the 
differences in the Lya line profiles predicted by the two outflow geometries in terms of their 
width, asymmetry and shift from the line centre for a set of outflows with different hydrogen 
column densities, expansion velocities and metallicities. In general, the Lya line profile of the 
Shell geometry is broader and more asymmetric, and the Lya escape fraction is lower than 
with the Wind geometry for the same set of parameters. In order to implement the outflow ge- 
ometries in the semi-analytical model GALFORM, a number of free parameters in the outflow 
model are set by matching the luminosity function of Lya emitters over the whole observed 
redshift range. The resulting neutral hydrogen column densities of the outflows for observed 
Lya emitters are predicted to be in the range ^ 10^^ — 10^'^[cm^^]. The models are con- 
sistent with the observationally inferred Lya escape fractions, equivalent width distributions 
and with the shape of the Lya line from composite spectra. Interestingly, our predicted UV 
luminosity function of Lya emitters and the fraction of Lya emitters in Lyman-break galaxy 
samples at high redshift are in partial agreement with observations. Attenuation of the Lya 
line by the presence of a neutral intergalactic medium at high redshift could be responsible 
for this disagreement. We predict that Lya emitters constitute a subset of the galaxy popu- 
lation with lower metallicities, lower instantaneous star formation rates and larger sizes than 
the overall population at the same UV luminosity. 
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1 INTRODUCTION 

Over the past 10 years the Lya line has proved to be a suc- 
cessful tracer of galaxies i n the redshift range 2 < z < 



7 (e. g . ICowie&H 



Kudritzki et ai] I2OOOI; iRhoads et alj 



200d: IHu et alj 12002 : Gronwall e t alj l2007i: lOuchi et alj 



200s 



20061 : 



Nilsso n et alj|2009l:'Shimasaku et alj |2006l " iKashikawa et alJ I 
Huet al. 2010; Guaita et al. 2010). More recently, samples of 

Lyg emitters at z ~ 0.2 obtained with the GA LEX satellite 

jPeharveng et alj|2008l : ICowie, Barger & Hij|2010l) . have allowed 
us to study this galaxy population over an even broader range of 
redshifts. 

Star-forming galaxies emit Lya radiation when ionizing pho- 
tons produced by massive young stars are absorbed by atomic hy- 
drogen (HI) regions in the interstellar medium (ISM). These hy- 
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drogen atoms then recombine leading to the emission of Lya pho- 
tons. Therefore, Lya emission is, in principle, closely related to the 
star formation rate (SFR) of galaxies. However, in general only a 
small fraction of L yg photons manage to escape from galaxies (e.g. 
iHaves et"al]|201lh . This makes it difficult to relate Lya emitters 
to other star forming galaxy populations at high redshift, such as 
Lyman-break galaxies (LBGs) or sub-millimetre galaxies (SMGs). 

The physical properties of galaxies selected by their 
Lya emission are inferred from spectral a nd photometric data 
jGawiser et alj l2007l : iGronwall et al] I200I iNilsson et all l2009l : 
I Guaita etaL 2011 ). Furthermore, Lya emitters are currently used 
to study the kinematics of the ISM i n hig h redshift galaxies 
jShaplev et al]|2003l : ISteidel et alj|2010l, I2OI1I : ^ulas et al j|201lh, 
to trace the large scale struc ture of the Universe UShimasaku et alj 

1] l2007l : lOrsi et alj l200a 
2010), to constrain the epoch 
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Despite this progress, understanding the physical mechanisms 
which drive the escape of Lya radiation from a galaxy remains a 
challenge. Lya photons undergo resonant scattering when interact- 
ing with hydrogen atoms, resulting in an increase of the path length 
that photons need to travel before escaping the medium. Therefore, 
the probability of photons being absorbed by dust grains is greatly 
enhanced, making the escape of Lya photons very sensitive to even 
small amounts of dust. Furthermore, as a result of the complex ra- 
diative transfer, the frequency of the escaping Lya photons gen- 
erally departs from the line centre as a consequence of the large 
number of scatterings with many hydrogen atoms. 

Recent observational studies have use d var ious methods to in- 
fer the escape frac ti on of Lya photons, fescjAtek et al. 2008, 2009|; 
lOstUn et al.l l2009l : iKornei et al.l I2OI0I : iHaves et alj I20ia I20lih . 
This is generally done either by comparing the observed line ratio 
between Lya and other hydrogen recombination lines, such as Ha 
and H/?, or by comparing the star formation rate derived from the 
Lya luminosity to that obtained from the ultraviolet continuum. 
The first method is the more direct, since the intrinsic fluxes of 
the comparison lines can be inferred after correcting the observed 
fluxes for extinction. Then, the departure from case B recombina- 
tion of the ratio of the Lya intensity to another hydrogen recom- 
bination line is attributed to the escape fraction of Lya differing 
from unity. The second method, on the other hand, relies heavily 
on the assumed stellar population model used, the choice of the 
stellar initial mass function (IMF) and the attenuation of the ultra- 
violet continuum by dust, and is therefore more uncertain. 

These measurements have revealed that the escape fraction of 
Lya emitters can be anything between 10~^ and 1. The observa- 
tional data listed above also suggest a correlation between the value 
of the escape fraction and the dust extinction. The large scatter 
found in this relation suggests there is a range of physical parame- 
ters which determine the value of feac • 

Early theoretical models of L ya emission from galaxies were 
based on a static ISM (see, e.g. iNeufeld [T990l : ICharlot & Fall 
1 19931) . These models explained the difficulty of observing Lya 
in emission due to its very high sensitivity to dust in such a 
medium. Moreover, the first observations of Lya emission in lo- 
cal starburst dwarf galaxies suggested a strong correlation be- 
tween metallicity and Lya luminosity 1 Meier & Terlevich|[l98ll : 
iHartmann. Huchra & Geile3ll984l : iHartmann et alJll988h . leading 
to the conclusion that metallicity, which supposedly traces the 
amount of dust in galaxies, is the most important factor driving the 
visibility of the Lya line. 

However, subsequent observational studies showed only a 
weak correlation between Lya luminosity and metallicity, sug- 
gesting instead the importance of the neutral gas distribution and 
its kinematics (e.g., Giavali sco. Koratkar & Calzettilll99 6). Fur- 
ther analysis of metal lines in local starbursts revealed the pres- 
ence of outflows which allow the escape of Lya photons. The ob- 
served asymmetric P-Cygni Lya line profiles are consistent with 
Lya photons escaping from an expand i ng shell of neutral ga s 
jThuMi & Izotovlll997flKunth et alJI 19981 ; iMas-Hesse e7ani2003h . 
This established outflows as the main mechanism responsible for 
the escape of Lya photons from galaxies. Furthermore, observa- 
tions at higher redshifts reveal Lya line profiles which also re- 
semble those expected when photons escape through a galactic 



outflow iShaplev et alj |2003| ; iKashikawa et alj I2OO6I : iKomei et al] 
l2010l : lHuetalJl2010h . 

In the last few years there has been significant progress in the 
modelling of Lya emitters in a cosmological setting. The first con- 
sistent hierarchical galaxy formati on model which i ncluded Lya 
emiss ion is th e one described by iLe Delliou et al. I ( I2OO5I I2OO6I) 
and Orsi et alj ( |2008|) . which makes use of the GALFORM semi- 
analytical model. In this model, the simple assumption of a fixed 
escape fraction, fesc = 0.02, regardless of any galaxy property or 
redshift, allowed us to predict remarkably well the abundances and 
clustering of Lya emitters over a wide range of redshifts and lumi- 
nosi ties. 

iNagamine et alj |2003.|2010|) modelled Lya emitters in cos- 
mological SPH simulations. In order to match the abundances 
of Lya emitters at different redshifts, they were forced to in- 
troduce a tunable escape fractio n and a dut y cycle parameter. 
[Kobavashi. Totani & Nagashimdl ( I200I I2OI0I) developed a sim- 
ple phenomenological model to compute fesc in a semianalyti- 
cal model. Their analytical prescription for fesc distinguishes be- 
tween outflows produced in starbursts and static media in qui- 
escent galaxies. Dayal et. al use an SPH simulation to study 
Lya emitters at high redshift and their attenuation by the neu- 
tral intergalactic medium (IGM). However, they assume the 
Lya esca pe fraction is related to the escape of UV continuum 
photonsJDayal, Ferrara & Gallerani 2008; Dayal, Ferrara & Sard 
l20ld ; [Pava l, Maselli & Ferrara 2011). Tilvi et al. (2009, 201^ 
make predictions for Lya emitters using an N-body simulation and 
the uncertain assumption that the Lya lumin osity is proportional to 
the ha lo mass accretion rate. More recently. iForero-Romero et aU 

ij) presented a model for high redshift Lya emitters based on 
a hydrodynamic simulation, approximating Lya photons to escape 
from homogeneous and clumpy gaseous, static slabs. 

Motivated by observational evidence showing that Lya pho- 
tons escape through outflows, we present a model that incorporates 
a more physical treatment of the Lya propagation than previous 
work, whilst at the same time being computationally efficient, so 
as to allow its application to a large sample of galaxies at different 
redshifts. 

Such a physical approach to modelling the escape of Lya 
photons requires a treatment of the radiative transfer of photons 
through an HI region. The scattering and destruction of Lya pho- 
tons have been extensi vely studied due to their many applications 
in astrophysical media. 'Harrin gtoril ( I1973I) studied analytically the 
emergent spectrum from an optically thick, homogeneous static 
slab with p hotons generated at the line centre. This result was gen- 
eralised by iNeufeld ( Il990l) . to include photons generated at any 
frequency, and to provide an analytical expression for fesc in this 
configuration. 

Numerical methods, on the other hand, allow us to study 
the line profiles and escape fractions of Lya photons in a 
wider variety of configurations. The standard approach is to 
use a Monte Carlo algorithm, in which the paths of a set of 
photons are followed one at a time through many scattering 
events, until the photon either escapes or is absorbed by a 
dust grain. Such calculations have been applied successfully 
to study the properties of Ly a emitters in different scenar- 

2000; Zheng & Miralda-Escudj 



ios ( see, 
2002 



(s ee, e . g. lAhn. Lee & Le ^ ,2000; Zheng & Miralda- Escuc 
ij; lAhnl I200I |2004 IVerhamme Schaerer & Masellil (200i 



Dii kstra. Haiman & Spaans 200d; Laursen & Sommer-LarsenI 



2007i; iLaursen.' Razoumov & Sommer-Larsenll2009l) ~ 

Recently, IZheng et alj feOld 1201 ih combined a Monte Carlo 
Lya radiative transfer model with a cosmological reionization sim- 
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ulation at 2: ~ 6, obtaining extended Lya emission due to spatial 
diffusion. Their simulation box is, however, too small to be evolved 
to the present day without density fluctuations becoming nonlinear 
on the box scale, and does not have a volume large enough to sam- 
ple a wide range of environments. 

Previous work has not studied Lya emitters in a framework 
that at the same time spans the galaxy formation and evolution pro- 
cess over a broad range of redshifts and includes Lya radiative 
transfer. The need for such hybrid approach is precisely the moti- 
vation for this paper. 

Given the observational evidence that outflows facilitate the 
escape of Lya photons from galaxies, here we study the nature 
of Lya emitters by computing the escape of photons from galax- 
ies in an outflow of material by using a Monte Carlo Lya radia- 
tive transfer model. Galactic outflows in our model are defined ac- 
cording to predicted galaxy properties in a simple way. This makes 
our modelling feasible on a cosmological scale, whilst retaining all 
the complexity of Lya radiative transfer. Following our previous 
work, we use the semianalytical model GALFORM. This paper rep- 
resents a significant improvement over the treatment of Lya emit- 
ters in hierarchical galaxy for mati on models initial ly described in 
iLe Delliou eFZI ( l2005ll2006h and lOrsi et ai] ( |2008|) . which all as- 
sumed a constant escape fraction. 

The outline of this paper is as follows: Section[2]describes the 
galaxy formation and radiative transfer models used. Also, we de- 
scribe the two versions of galactic outflow geometries we use and 
explain how these are constructed in terms of galaxy parameters 
that can be extracted from our semi-analytical model. In Section[3] 
we present the properties of the Lya line profiles and escape frac- 
tions predicted by our outflow geometries. In Section|4]we describe 
the galaxy properties that are relevant to our Lya modelling and 
describe how we combine the Lya radiative transfer model with 
the semi-analytical model. In Section [5] we present our main re- 
sults, where we compare with observational measurements when- 
ever possible, spanning the redshift range < 2 < 7. In Sec- 
tion |6] we present the implications of our modelling for the prop- 
erties of galaxies selected by their Lya emission, compared to the 
bulk of the galaxy population. Finally, in SectionQwe discuss our 
main findings. In the Appendix we give details of the implementa- 
tion of the Monte Carlo radiative transfer model, the effect of the 
UV background on the outflow geometries and the numerical strat- 
egy followed to compute the escape fraction of Lya photons from 
galaxies. 



2 MODEL DESCRIPTION 

Our approach to modelling Lya emitters in a cosmological frame- 
work involves combining two independent codes. The backbone of 
our calculation is the GALFORM semi-analytical model of galaxy 
formation, outlined in Section [ZT] from which all relevant galaxy 
properties can be obtained, including the intrinsic Lya luminosity. 
The second is the Monte Carlo Lya radiative transfer code to com- 
pute both the frequency distribution of Lya photons and the Lya 
escape fraction. This code is described briefly in Section |Z21 with 
further information and tests presented in Appendix lAl In Section 
|2.3| we describe the outflow geometries and the way these are de- 
fined in terms of galaxy properties. 



2.1 Galaxy formation model 

We use the semi-analytical model of galaxy formation, GALFORM, 
to predict the properties of galaxies and their abundance as a 
function of redshift. The GALFORM mod el is fully described in 
Cole et al. (2000) and lBenson et alj j2003h (see also t he review by 
lBaugh.2006) . The variant used here was i ntroduced bylSaugh et alj 
( l2005h . and is also described in detail in iLacev et alj j2008l) . The 
model computes star formation and galaxy merger histories for the 
whole galaxy population, following the hierarchical evolution of 
the host dark matter haloes. 

The iBaugh et all feOOSh model used here is the same 
GALFORM variant used in our previous wor k on Lya emitters 
jLe Delliou e"tar "2005. 2006: Orsi et alj |2008h . A critical assump- 
tion of the iBaugh et al. (2005) model is that stars which form 
in starbursts have a top-heavy initial mass function (IMF). The 
IMF is given by dA'^/dln(m) cx: with x = in this case. 

Stars formed quiescently in discs have a solar neighbourhood IMF, 
with the form proposed by iKennicutJ l ll983l) with x = 0.4 for 
m < IA/0 and x = 1.5 for m > IMq. Both IMFs cover the 
mass range O.15M0 < m < 125M0. Within the framework of 
ACDM, Baugh et al. argued that the top-heavy IMF is essential 
to match the counts and redshift distribution of galaxies detected 
through their sub-millimetre emission, whilst retaining the match 
to galaxy properties in the local Universe, such as the optical and 
far-IR luminosity f unctions, galaxy gas fractions and metallicities. 
iLacev et"ail ^2008 ) showed that the model predicts galaxy evolu- 
tion in the IR in good agreem ent with observations from Spitzer. 
Moreover, iLacev et au j201lh also showed that the Baugh et al. 
model predicts the abundance of Lyman-break galaxies (LBGs) 
remarkably well ove r the redshift range 3 < 2 < 10 (see also 
lGonzalezetalj20Ilh . 

In GALFORM, the suppression of gas cooling from ionizing 
radiation produced by stars and active galactic nuclei (AGNs) dur- 
ing the epoch of reionization is modelled in a simple way: after the 
redshift of reionization, taken to be Zrdon = 10, photoionization of 
the IGM completely suppresses the cooling and collapse of gas in 
haloes with circular velocities below Krit ~ 30[km/s]. 

The Baugh et al. model assumes two distinct modes of feed- 
back from supernovae, a reheating mode, in which cold gas is ex- 
pelled back to the hot halo, and a superwind mode, in which cold 
gas is ejected out of the galaxy halo. We describe both modes 
of supernova fee dback in more detail in section 12.3.21 (see also 
iLacev et al.ll2008h . 

Unlike the iBower et J] ( I2OO6I) variant of GALFORM, the 
model used here does not incorporate feedback from an AGN. 
The superwind mode of feedback produces similar consequences 
to the quenching of gas cooling by the action of an AGN, as 
both mechanisms s uppress the bright end of the LF. However, the 
iBower et al.l ( l2006h model does not reproduce the observed abun- 
dances of LBGs or submillimetre galaxies (SMGs) at high redshift, 
which span a redshift range which overlaps with that of the LFs 
of Lya emitters considered in this paper. Therefore, we do not use 
this GALFORM variant to make predictions for Lya emitters. 

In GALFORM, the intrinsic Lya emission is computed as fol- 
lows: 

(i) The composite stellar spectrum of the galaxy is calculated, 
based on its predicted star formation history, including the effect 
of the metallicity with which new stars are formed, and taking into 
account the IMFs adopted for different modes of star formation. 

(ii) The rate of production of Lyman continuum (Lyc) photons 
is computed by integrating over the composite stellar spectrum. We 
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assume that all of these ionizing photons are absorbed by neutral 
hydrogen within the galaxy. The mean number of Lya photons pro- 
duced per Lyc photon, assuming a gas temperature of 10* K for the 
ioniz ed gas, is approxi mately 0.677, assuming Case B recombina- 
tion dOsterbrockll 19891) . 

(iii) The observed Lya flux depends on the fraction of Lya 
photons which escape from the galaxy (/esc). Previously we as- 
sumed this to be constant and independent of galaxy properties. 
The escape fraction was fixed at fcsc = 0.02, resulting in a remark- 
ably good match to the observed Lya LFs over the redshift range 
3 < 2 < 6. In our new model, we make use of a Monte Carlo 
radiative transfer model of Lya photons to compute fesc in a more 
physical way. This model is briefly described in the next subsec- 
tion. 



2.2 Monte Carlo Radiative Transfer Model 

In order to compute the Lya escape fraction and line profile we 
construct a Monte Carlo radiative transfer model for Lya photons. 
This simulates the escape of a set of photons from a source as they 
travel through an expanding HI region, which may contain dust, by 
following the scattering histories of individual photons. By follow- 
ing a large number of photons we can compute properties such as 
the Lya spectrum and the escape fraction. 

Our Monte Carlo radiative transfer code works on a 3D grid 
in which each cell contains information about the neutral hydrogen 
density, uh, the temperature of the gas, T, and the bulk velocity, 
fbuik. Once a Lya photon is created, a random direction and fre- 
quency are assigned to it, and the code follows its trajectory and 
computes each scattering event of the photon as it crosses the HI 
region until it either escapes or is absorbed by a dust grain. If the 
photon escapes, then its final frequency is recorded, a new photon 
is generated and the procedure is repeated. 

Clearly, to get an accurate description of the escape of Lya 
photons from a given geometrical setup, many photons must be fol- 
lowed. The number of photons needed to achieve convergence will 
depend on the properties of the medium, but also on the quantity 
in which we are interested. For example, for most of the outflows 
studied here, only a few thousand photons are needed to compute 
an accurate escape fraction. However, tens of thousand of photons 
are needed to obtain a smooth line profile. 

Our radiative transfer model of Lya photons 
is similar to previous models in the literature (e. 

Ahn 



Zheng & Miralda-Escude 



2002; 



2003', 20041: 



Verhamme, Schaerer & Maselli 2006; Dijkstra, Haiman & SpaansI 
20061; IL aursen, Ra zoumov & Sommer-LarsenI 1 20091 : 

Barnes & Haehnelt .20101) . We describe the numerical imple- 



mentation of our Monte Carlo model and the validation tests 
applied to it in Appendix lAl 



2.3 Outflow geometries 

In our model, the physical conditions in the medium used to com- 
pute the escape of Lya photons depend on several properties of 
galaxies predicted by GALFORM. Below we describe two outflow 
geometries for the HI region surrounding the source of Lya pho- 
tons. They differ in their geometry and the way the properties of 
galaxies from GALFORM are used. We assume the temperature of 
the medium to be constant at T = 10* K, which sets the thermal 
velocity dispersion of atoms, following a Maxwell-Boltzmann dis- 
tribution, to be vt\i = 12.85 kms~^ (see Eq. lAll in Appendix lAlfor 



more details). For simplicity, the source of Lya photons is taken to 
be at rest in the frame of the galaxy, in the centre of the outflow, 
and emits photons at the line centre only, A = 1215.68 A. 

2.3.1 Expanding thin shell 

Previous radiative transfer studies of Lya line profiles adopted 
an expanding shell in the same way as we used here (see, e.g. 
Ahn 2003, 2004; Verham me. Schaerer & Masellill2006l : ISchaereJ 
|2007| ; iVerha mme et al.i i20()^ and references therein). This model, 
hereafter named "Shell", consists of a homogeneous, expanding, 
isothermal spherical shell, in which dust and gas are uniformly 
mixed. The shell, although thin, is described by an inner and 
outer radius i?in and i?out, which satisfy i?in = /thi?out. We set 
/th = 0.9. In addition, the medium is assumed to be expanding 
radially with constant velocity Vixp. The column density through 
the Shell is given by 



NH{r) 



(1) 



ATxruHRl^t ' 

where Xh ~ 0.74 is the fraction of hydrogen in the cold gas and 
tuh is the mass of a hydrogen atom. In GALFORM, the Lya lumi- 
nosity originates in the disk (in quiescent galaxies) or the bulge (in 
starbursts). Some galaxies may also have contributions from both 
components. Therefore, A/shdi, Rout and Kxp are taken to be pro- 
portional to the mass of cold gas AfcoW, half-mass radius R1/2 and 
circular velocity Vcirc, respectively, i.e.. 



-Rout = /fl{-Rl/2), 



14> 



(2) 
(3) 
(4) 



where /a/, /r and fv are free parameters. To take into account the 
contribution from both components of a galaxy, we define 

(A/gas) = i^Lt^A^gif + (1 - -FLi^)A^gas'"^ (5) 

I r> \ _ pdisk pdisk I pdisks pbulge ^^^^ 

\-tt-l/2/ — -f^Lya-rtl/2 + ^ ^ Lya ) /2 i ("> 

(7) 



{vii.c) = i^Lt^K?;? + (1 - FLt')Ki"c''^ 



j-idisk 
^Lya 



T disk 
T total 



(8) 



In each term, the superscript indicates the contribution from the 
disk, the bulge or the total (the sum of the two). In most galaxies, 
however, either the disk or bulge term completely dominates. 

Likewise, the metallicity of the shell Zout is taken to be the 
metallicity of the cold gas ZcoM weighted by a combination of the 
mass of cold gas and the Lya luminosity of each component, i.e., 

Mdisk IT disk ^disk , ii ^bulge r bulge ybulge 
gas -i^Lya^cold + -"'Jgas ^Lya '^cold 



Zout 



/\,f disk r disk I n^bulgo bulge 
JWgas ^hya + -'"gas ^hya 



(9) 



In order to compute the dust content in the outflow we assume that 
the mass of dust in the outflow. A/dust, is proportional to the gas 
mass and metallicity, i.e. 



(5. 

Afdust ~ 75— Afshcll^out, 

Zq 



(10) 



where the dust-to-gas r atio is taken to be S t = 1/110 at the solar 
metallicity Zq = 0.02 dOranato et alj2000l) . The extinction optical 
depth at the wavelength of Lya can be written as 



Td 



-^^NnZont, 

Zq 



(11) 
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where Eq = 1.77 x 10"^^ [cm^] is the ratio ra/Nn for so- 
lar metallicity at the wavelength of Lya (1216 A). Throughout 
this w ork we use the extinction curve and albedo from Silva et al. 
( Il998l) . which are fit to the observed extinction and albedo in the 
Galactic ISM. For a dust albedo A, the optical depth for absorption 
is simply 

= (1 - A)Td. (12) 
At the wavelength of Lya, ALya = 0.39. 

2.3.2 Galactic Wind 

Supemovae heat and accelerate the ISM th rough shock s and hence 
gener ate outflows from galaxies (see, e.g. lFrank|[l99^ : IStricklan j 
I2OO2I) . Here we develop an outflow model, hereafter called "Wind", 
which relates the density of the outflow to the mass ejection rate 
from galaxies due to supemovae. In GALFORM, this mass ejection 
rate is given by 



3roh(K;irc) + ;flsw(Vcirc)] V": 



where 

Pich 
/3sw 



T/ \ -a hot 

/swmin[l, (Kirc/14w)' 



(13) 

(14) 
(15) 



The terms /3rch and /3sw define the two different modes of super- 
nova feedback (the reheating and superwind), and the constants 
Qhot, Vhot, /sjD and Vsw are free parameters of GALFORM, cho- 
sen by fitting the model predictions to observed galaxy LPs. The 
instantaneous star formation rate, ■(/) is obtained as 



A4, 



(16) 



where r* is the star formation time-scale, which is different for 
quiescent galaxies and starbursts. For a detailed description of the 
star formation and supernova feedb ack processes in t his va riant of 
GALFORM, see lBaugh et alj |2005') and Lacev et alj ( I2OO8D . Since 
star formation can occur in the disk and in the bulge, there is a mass 
ejection rate term Mej associated with each component. 

We construct the wind as an isothermal, spherical flow ex- 
panding at constant velocity Kxp, and inner radius -Rwind (the wind 
is empty inside). In a steady-state spherical wind, the mass ejection 
rate is related to the velocity and density at any point of the Wind 
via the equation of mass continuity, i.e. 



Mej = 47rrVcxpp(r) 



(17) 



where p(r) is the mass density of the medium, and Vcxp is calcu- 
lated following equation l|4j. Since star formation in GALFORM can 
occur in both the disk and bulge of galaxies, A/ej in equation (TTJl 
corresponds to the sum of the ejection rate from the disk and the 
bulge. 

Following equation l |17| l, the number density profile nHi{r) 
in the Wind geometry varies according to 







riHiir) 



The column density Nh of the outflow is 



r < R 

Wind 
r ^ R-wind- 



Nh 



47rmH-Rwind Vcxp 



(18) 



(19) 



where the inner radius of the wind, i?wind, is computed in an anal- 
ogous way to -Rout in the Shell geometry (Eq.[3j. Note that both fn 
and fv in this case are different free parameters and independent 
of those used in the Shell geometry. The metallicity of the Wind, 
Zviind, corresponds to the metallicity of the cold gas in the disk and 
bulge weighted by their respective mass ejection rates. 

For computational reasons, the radiative transfer code requires 
us to define an outer radius, -Router, for the Wind. However, since 



the number density of atoms decreases as 



we expect that 



at some point away from the centre the optical depth becomes so 
small that the photons will be able to escape regardless of the ex- 
act extent of the outflow. We find that an outer radius -Router ~ 
20-Rwind yields converged results, i.e. the escape fraction of Lya 
photons does not vary if we increase the value of -Router further. 



3 OUTFLOW PROPERTIES 

In this section we explore the properties of the Lya radiative trans- 
fer in our outflow geometries prior to coupling these to the galaxy 
formation model GALFORM. We do this by running our Monte 
Carlo radiative transfer model over a grid of configurations span- 
ning a wide range of hydrogen column densities, expansion veloci- 
ties and metallicities. In order to obtain well-defined Lya line pro- 
files, we run each configuration using 5 x 10'* photons. Therefore, 
the minimum Lya escape fraction our models can compute in this 
case is fcse = 2 x 10^^. 



3.1 Lya line properties 

One important difference between our two outflow geometries is 
the way the column densities are computed. In the Shell geometry 
the column density is a function of the total cold gas mass of a 
galaxy and its size, whereas in the Wind geometry it depends on 
the mass ejection rate given by the supernova feedback model along 
with the size and the circular velocity of the galaxy. This difference 
translates into different predicted properties for Lya emitters, as 
will be shown in the next section. However, even when the column 
density, velocity of expansion and outflow metallicity are the same, 
the two models will give different escape fractions and line profile 
shapes due to their different geometries. 

Fig.[T]shows the Lya line profiles obtained with the two mod- 
els when matching the key properties for outflows at the same col- 
umn densities. In order to make a fair comparison between the Shell 
and Wind geometries, we compare configurations with the same 
column density, expansion velocity and metallicity. In addition, the 
inner radius in the Wind geometry is chosen to be the same as its 
counterpart, the outer radius in the Shell geometry. The two pan- 
els of Fig. [T] display a set of configurations with column densities 
of Nh = 10^" [cm-^] (left panel) and iV^ = 10^^ [cm"^] (right 
panel). In the following, we express the photon's frequency in terms 
of X, the shift in frequency around the line centre uo, in units of the 
thermal width, i.e. 



(20) 



where Ai^n = uthJ^o/c, and c is the speed of light, and uth is the 
thermal velocity discussed in Section [231 

As a general result, outflows with a column density of Nh ~ 
lO^^lcm"^], regardless of the other properties, produce multiple 
peaks at frequencies redward of the line centre in both geometries, 
with different levels of asymmetry. The Shell geometry generates 
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Figure 1. Comparison of the Lya line profile obtained witli the Wind and Shell geometries for Nh = lO'^" [cm~'^] (leftmost two columns) and Nj] = 
lO'^^ [cm~'^] (rightmost two columns). The red (blue) histogram shows the full spectrum obtained for the Wind (Shell) outflow geometry. The cyan, green, 
pink and gray lines show the spectra of photons which experienced 0,1,2 and 3 or more back-scatterings before escaping, respectively. Likewise, the same 
colours are used to plot vertical dashed lines showing the frequencies corresponding to —xi,s, —2x)jg, —3xbs and ~4xbg. respectively (see text for details). 
Each row displays a different configuration, characterised by a given Nu, Vcxp and Z, as indicated in the legend. The top two rows con'espond to dust-free 
configurations (Z = 0), whereas the bottom two have different metallicities, chosen to have equal dust optical depth, although having different column 
densities of hydrogen. The escape fraction of Lya photons, fcsc, is indicated in each box. The Lya profiles shown are normalized to the total number of 
escaping photons (instead of the total number of photons run), to ease the comparison between the dusty and dust-free cases. Note the different range in the 
X-axis between the left and right panels. 



more prominent peaks than the Wind geometry. The frequency of 
the main peak is, however, the same in both outflow geometries. 

On the other hand, outflows with Nh = 10^^[cm~^] display 
broader Lya profiles. As opposed to the lower column density case, 
these profiles display a single peak, also redward of the line centre. 
The position of this peak is also the same for both geometries. 

The effect of a large expansion velocity is shown in the 2nd 
and 4th rows in both panels of Fig.[T] Here the configurations have 
Vexp = 500[km/s]. The differences with the configurations where 
T4xp ~ 200[km/s] are evident: the profiles are broader, and the 
position of the peaks are displaced to redder frequencies. Further- 
more, in the left panel of Fig.[I] where Nh = 10^"[cni~^], it is 



shown that a significant fraction of photons escape at the line cen- 
tre. The high expansion velocity in this case makes the medium op- 
tically thin, allowing many photons to escape without undergoing 
any scattering. 

For the configurations with Nh ~ lO^^fcm"^], the optical 
depth that a photon at the line centre (x = 0) sees when traveling 
along the radial direction is tq — 3.63 and to = 0.57 for expan- 
sion velocities Vaxp ~ 200[km/s] and T4xp = 500[km/s], respec- 
tively. Accordingly, the configurations with Nh = 10^'^[cm~^] 
have optical depths a factor 100 higher. On the other hand, a static 
medium with Nh = 10^" [cm~^] has a much higher optical depth, 
To = 3.31 X lO*". This illustrates the strong effect of the expansion 
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Figure 2. The properties of the Lyo profiles predicted by the two outflow geometries for a given column density, expansion velocity and metallicity. Left 
panels: the width of the profiles, defined as the difference between the 90 and 10 percentiles of the frequency distribution of escaping Lya photons, measured 
in km/s. Middle panels: the asymmetry, defined in the text, and measured in km/s. Right panels: the shift of the median of the Ly« profiles with respect to 
the line centre, measured in km/s. Note that negative velocities indicate redshifting. The top row shows the Lyo profile propeilies in dust-free outflows. The 
bottom row shows outflows with a metallicity Z = 0.01. Different expansion velocities spanning the range 100 — 500 km/s are shown with different colours, 
as shown by the labels in each panel. The Wind geometry is shown with solid lines and filled circles. The Shell geometry is shown with dashed lines and open 
circles. 



velocity in reducing the optical depth of the medium, thus allowing 
photons to escape. 

The Lya line profiles obtained can be characterised by the 
frequency distribution of photons split according to the number of 
backscatterings they experience before escaping (i.e. the number of 
times photons bounce back to the inner, empty region). When pho- 
tons interact for the first time with the outflow, a fraction of them 
will experience a backscattering. These events are significant since 
the distribution of scattering angles is dipolar (see Eq. |A13| in Ap- 
pendix [A]). The frequency of a photon after a scattering event, in 
the observer's frame, is given by Eq. dAlSb . Depending on the di- 
rection of the photon after the scattering event, its frequency will 
fall within the range x — \—2xi,s,0], where Xbs = Vexp/t>th (see 
also lAhnl2063l : IVerhamme, Schaerer & Masellil26o3) . Photons that 
do not experience a backscattering, or escape directly, form the 
cyan curves in Fig.[T] If the photon is backscattered exactly back- 
wards, then its frequency will be a; = —2xbs- The cross-section 
for scattering is significantly reduced for these photons, and so a 
fraction escape without undergoing any further interaction with the 
outflow. Photons backscattered once form the green curves in Fig. 
[T] If a photon experiences a second backscatter in the exact oppo- 
site direction, then its frequency will become x — —Axts, and the 
cross-section for a further scattering will again reduce significantly. 
The magenta curves show the distribution of photons that experi- 
enced 2 backscatterings. Finally, photons that experience 3 or more 
backscatterings are shown in gray. 



In detail, the geometrical differences between the Shell and 
Wind models are translated into each backscattering peak con- 
tributing in a different proportion and with a different shape to 
the overall spectrum for each geometry. Previous studies have also 
found this relation between the peaks of backscattered photons 
and a;bs in m edia with c o lumn d ensities of the order of Nh ^ 
jAhi] |2003| , |2004 IVerhamme, Schaerer & Masellil 
l2006t) , although they did not study the line profiles for higher op- 
tical depths as we do here. For Nh = 10^'^[cm~^], we find that 
the peaks are displaced considerably from their expected position 
based on the simple argument above. This is not surprising, since 
in outflows with very large optical depths the number of scatterings 
broadens the profiles and reddens the peak positions. 

In the Wind geometry, the contribution of photons featuring 
no backscatterings dominates most of the total profile for Nh = 
10'^" [cm~^], whereas in the overall line profile for the Shell ge- 
ometry there is a clear distinction between a region dominated 
by photons suffering no backscatterings and those backscattered 
once (green curve). This illustrates again how the Lya line pro- 
file in the Shell geometry features clear multiple peaks from one or 
more backscatterings, whereas in the Wind geometry the secondary 
peaks are less obvious. 

Fig.[T]also shows the effect of including dust in the outflows. 
Overall, dust absorption has more effect on the redder side of the 
line profiles than at frequencies closer to the line centre, where 
the probability of scattering with hydrogen atoms is significantly 
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higher than the probability of interaction with dust. In detail, the 
sensitivity of Lya photons to dust is the result of an interplay 
between the optical depths of hydrogen and dust at the photon 
frequency (the cross-section of scattering is significantly reduced 
away from the line centre, and hence so is the average number of 
scatterings) and the bulk velocity of the gas, reducing the number 
of scatterings further, and thus also the probability for a photon to 
be absorbed. 

A less obvious outflow property for determining the Lya at- 
tenuation by dust is its geometry. Despite the similarities between 
the outflow geometries, the Shell geometry consistently gives a 
lower Lya escape fraction than the Wind geometry in the configu- 
rations studied in Fig.[T] 

As a final comparison, in Fig. [T] we chose a metallicity of 
Z = 0.2 for the configurations with Nh ~ lO^^fcm"^] and 
Z = 0.002 for the configurations with Nh = 10^^ [cm"^]. Al- 
though the metallicities are different, the optical depth of dust, 
given by Eq. i ll It . is the same in both cases, Td = 1.77. By match- 
ing the optical depth of dust, we can study the effect of Nh on 
the Lya escape fraction. Fig. [T] shows that even when the optical 
depth of dust is the same, outflows with A'^^r = 10^^[cm~^] have 
Lya escape fractions up to an order of magnitude lower than those 
with Nh = 10'^''[cm~'^]. This occurs since, in the former case, 
the average number of scatterings is about two orders of magnitude 
larger than in the latter, and hence the probability of photons inter- 
acting with a dust grain increases accordingly. This illustrates the 
complexity of the Lya radiative transfer process. 

To gain more insight into the difference in the Lya profiles 
generated by different configurations, we show in Fig. |2] a mea- 
sure of the width, asymmetry and average frequency shift of the 
Lya profiles for a set of configurations spanning a range of ex- 
pansion velocities of V'exp[km/s] = 100 — 500, column densities 
A^iffcm"^] = 10^^ - 10^^, and also two metallicities, Z = 
and Z = 0.01. To measure the width of the profiles we compute 
the difference in frequency between the 90 and 10 percentiles, ex- 
pressed in velocity units. The asymmetry is computed as the differ- 
ence between the blue side of the profile, Pgo — P50 and the red 
side, P50 — Pio, where Pio, P50 and P90 are the 10, 50 (median) 
and 90 percentile of the frequency distribution. The shift of the me- 
dian is simply the position of the median of the Lya profiles in 
velocity units. 

Overall, by examining the top and bottom rows of FiglHit be- 
comes clear that the width, asymmetry and median shift of the Lya 
profiles are fairly insensitive to the presence of dust in both outflow 
geometries. The only important effect of dust is to limit the range of 
column densities where the Lya profiles are appreciable. Configu- 
rations with log(7V£f [cm~'^]) > 21.5 do not feature data points on 
the plot, since all photons used to compute the Lya profiles were 
absorbed by dust in this case. The optical depth of absorption, Ta, 
increases in proportion with the column density, as shown by Eqs. 
(11) and 

The small contribution of dust in shaping the Lya profiles lim- 
its the information that can be extracted observationally from the 
Lya line profile. If the outflow geometries studied here are good 
approximations to the outflows in Lya emitters, then the Lya es- 
cape fraction cannot be usefully constrained by using the shape of 
the spectrum. 

Despite the above, the shape of the Lya profiles in our out- 
flow geometries is sensitive to other properties. Fig.[2]shows a clear 
increase in the width of the Lya profiles, with both increasing col- 
umn density and expansion velocity. Also, the Shell geometry gen- 
erates consistently broader profiles than the Wind geometry, for all 
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Figure 3. The Lya escape fraction predicted by the two outflow geome- 
tries as a function of column density, expansion velocity and metallicity. 
The top panel shows configurations with Z = 0.004. In the bottom panel 
Z = 0.02. Expansion velocities spanning the range 100 — 500 km/s ai'e 
shown with different colours, as shown in the legend of each box. The Wind 
geometry is shown with solid lines and filled circles. The Shell geometry is 
shown with dashed lines and open circles. The gray circles show the result 
of switching off the scattering of Lya photons due to hydrogen atoms in 
the model, making photons interact only with dust grains. 



the configurations studied. Note that this is illustrated in Fig. [T] as 
well. 

The broadening of the Lya profiles is, however, less evident 
in the configurations with Nh = 10^^ [cm^^]. At this low column 
density, and even at the lower expansion velocities, many photons 
manage to escape at the line centre without being scattered due 
to the reduced optical depth, as explained above. Hence, the Lya 
profiles are narrower. In particular, the width of the Lya profile, 
when Nh ~ 10^^[cm~^] and Kxp = 500[km/s] is zero, showing 
that nearly all photons escaped directly. 

At higher column densities, the expansion velocity plays an 
increasingly important role at broadening the profiles. In the out- 



flows with Nh = 10 
files increases from ~ 
velocities of Ve> 



[cm~ ], without dust, the width of the pro- 
300 [km/s] to ~ 1000 [km/s] for expansion 
100[km/s] and Voxp = 500[km/s], respec- 
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lively. When Nh = 10^^[cm~^], the width of the profile depends 
strongly on the geometry of the outflow. The width for the Wind 
geometry spans a range of ~ 1000[km/s] to ~ 2200[km/s] for ex- 
pansion velocities of Kxp = 100[km/s] and 14xp = 500[km/s], 
respectively, whereas for the Shell geometry it spans a much larger 
range of ~ 1700[km/s] to ~ 3200[km/s] for the same range of 
expansion velocities. 

The effect of dust on the width of the profiles is negligible, 
except for the Shell geometry with Nh = lO'^^ [cm~'^], where the 
width is reduced and is, therefore, closer to the width of the profiles 
of the Wind geometry. 

The asymmetry of the profiles is somewhat more complicated 
to characterise in terms of the column density and expansion veloc- 
ity of the outflows. As a general result, the Lya profiles are asym- 
metric towards the red-side of the spectrum. This is not surprising, 
since atoms in an outflow "see" the Lya photons redshifted, and the 
change in the direction of the photons due to the scattering event 
makes photons appear redder in the observer's frame. Overall, as 
Fig.|2]shows, the asymmetry is larger in the configurations with the 
highest expansion velocities when A'^H ^ 10^" [cm^^]. 

Finally, Fig. [2] also shows the shift of the median of the pro- 
files for the two outflow geometries. As expected, the median is 
redder with increasing column density and expansion velocity. The 
column density has a greater impact than the expansion velocity, 
since higher values of Nh imply a larger number of scatterings, 
thus increasing the reddening of the profiles. 



3.2 Escape fractions 

Fig.[3]compares the predicted Lya escape fractions in both outflow 
geometries for a set of column densities, expansion velocities and 
metallicities, similar to those considered in Fig. [2] 

As expected, both outflow geometries predict that the Lya es- 
cape fraction decreases rapidly with increasing Nh , as the medium 
becomes optically thicker. High values of the expansion velocity 
reduce the optical depth of the medium, hence enhancing the Lya 
escape fraction. 

For the range of properties studied here, we find that the Shell 
geometry predicts consistently lower Lya escape fractions than the 
Wind geometry, for the same set of parameters. This demonstrates 
that in the detailed interplay of physical conditions shaping fcsc, the 
outflow geometry plays an important role. 

The difference is less obvious in outflows with Nh < 
10^^ [cm^^], as is the influence of different expansion velocities. 
At larger column densities, even slightly different expansion veloc- 
ities can lead to significant differences in the resulting Lya escape 
fraction. 

Also, in Fig[3]we show the effect on fcsc of removing the scat- 
tering of photons by hydrogen atoms. We achieve this setting the 
Lya scattering cross section to zero as well (i.e. H(x) = in Eq. 
I A4I of Appendix [A^. The predicted fcsc is much higher than when 
considering the scattering by H atoms, and is also independent of 
expansion velocity, since in our modelling the optical depth of dust 
depends only on the metallicity and the column density of hydro- 
gen, as shown in Eq. Jllb . Fig[3] shows the effect of the resonant 
scattering resulting from the high cross-section at the line centre in- 
creasing the path length, and hence makes the resulting fcsc lower. 

It is interesting to note that even in the case of removing the 
scatterings by H atoms, the Shell geometry is more sensitive to dust 
than the Wind geometry, hence showing again the key role of the 
geometry of the outflows in determining the Lya fesc. 



4 THE MODEL FOR Lya EMITTERS 

In order to understand the nature of the predictions of our model 
for Lya emitters, we study in this section the galaxy properties pre- 
dicted by GALFORM that are relevant in determining the properties 
of Lya emitters (Section|4jTJ, prior to describing how we combine 
the radiative transfer model for Lya photons with GALFORM (Sec- 
tion|42t. 



4.1 Galaxy properties 

Each of our outflow geometries requires a series of galaxy prop- 
erties, provided by GALFORM, to compute the escape of Lya pho- 
tons. We consider the contribution of the disk and the bulge to com- 
pute averaged quantities, as described above. These are (i) the half- 
mass radius, R1/2, (ii) the circular velocity, Kirc, (iii) the metal- 
licity of the cold gas, ^coid, and (iv) the mass of cold gas of the 
galaxy, i\fgas. In addition, the Wind geometry requires the mass 
ejection rate due to supemovae, A/cj (see Eqns.[T3]and|17t. 

Fig.|4]shows the evolution of the galaxy properties listed above 
in the redshift range < z < 7, as a function of the intrinsic 
Lya luminosity, LLya,o- It is worth noting that these properties are 
extracted directly from GALFORM, so they do not depend on the 
details of the outflow model. 

In Fig. |4{a) we show the fraction of starbursts as a function 
of intrinsic Lya emission. Naturally, given the form of the IMF 
adopted, starbursts dominate in the brightest luminosity bins, re- 
gardless of redshift. However, the transition between quiescent and 
starburst Lya emitters shifts towards fainter luminosities as we go 
to lower redshifts. An important consequence of this trend is that 
the nature of Lya emitters, even at a fixed intrinsic Lya luminosity, 
is redshift dependent. Moreover, one might expect that the environ- 
ment in which Lya photons escape in quiescent galaxies and star- 
bursts is different. Although our fiducial outflow geometries do not 
make such a distinction, in section |42l we make our outflow geome- 
tries scale differently with redshift depending on whether galaxies 
are quiescent or starbursts. 

In Fig. |4fb) we show the dependence of the cold gas mass on 
the intrinsic Lya luminosity. As expected, in general the cold gas 
mass increases with LLya,o at a given redshift, since the latter is 
directly proportional to the star formation rate of galaxies, and, in 
this variant of GALFORM, the star formation rate is directly pro- 
portional to the cold gas mass, as shown in Eq. l ll6b . Note the star 
for mation timescale is different for quiescent and starburst galaxies 
(see lBaugh et alj|2005l for details on how the star formation is cal- 
culated in this variant of GALFORM, and iLagos et al.ll201ll for an 
alternative model). 

At low LLya,o quiescent galaxies dominate, whereas bright 
galaxies are predominately starbursts. However, there is a lumi- 
nosity range in which quiescent galaxies and starbursts both con- 
tribute. This is shown in Fig.|4tb) by a break in the cold gas mass- 
luminosity relation, which occurs in the luminosity range where 
low mass starbursts are as common as massive quiescent galax- 
ies. This luminosity corresponds to LLya[ergs~^ ~ 10*"^ 
at 2 ~ 0, and it shifts towards fainter Lya luminosities at higher 
redshifts. At 2 ~ 6, both quiescent and starbursts contribute at 

The total (disk + bulge) mass ejection rate is found to cor- 
relate strongly with the intrinsic Lya luminosity, as shown in Fig. 
|4{c). However, no significant evolution is found with redshift. Since 
the mass ejection rate due to supernovae is directly proportional to 
the star formation rate, as shown in Eqns. ( 114b and dlSt . and the 
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Figure 4. The evolution of the galaxy properties predicted by GALFORM which are used as inputs to the outflow geometries for z = 0.2 (green), z = 3.0 
(orange) and 2 = 6.6 (blue). The panels show, as a function of intrinsic Lyo luminosity, median values for (a) the fraction of starbursts; (b) the cold gas mass 
of the galaxy, M^oid; (c) the mass ejection rate, A/ej; (d) the luminosity-weighted half-mass radius, {R1/2)'' (s) 'he luminosity-weighted circular velocity, 
(Vcirc) ™d (f) the (mass and Lya luminosity)-weighted metallicity of the gas (.^gas) (Eq.|9). En'or bars show the 10-90 percentile range of the predicted 
distributions. 



conversion between the star formation rate and LLya,o depends on 
the production rate of Lyc photons, but not on redshift, it is not 
surprising that the mass ejection rate does not evolve with redshift. 

The half-mass radius {R1/2) is the parameter that is found to 
have the strongest evolution with redshift, as shown in Fig. IHd). 
Galaxies at z — 0.2 typically have half-mass radii of a few kpc/h. 
The median size of galaxies decreases rapidly with increasing red- 
shift, falling by an order of magnitude or more by 2 = 6.6. In 
contrast, {R1/2} varies only weakly with LLya.o- 

The circular velocity of galaxies depends both on their mass 
and half-mass radius. The sizes of galaxies are computed based 
on angular momentum conservation, centrigugal equilibrium for 
disks and virial equilibrium and en ergy conservation in mergers for 
spheroids, as described in detail bv lCole et alj j200Ct) . We find that 
the sizes correlate only weakly with the intrinsic Lya luminosity, 
as shown in Fig.|4td), and although the total mass of the galaxy is 
not necessarily related to the cold gas mass, the circular velocity 
(VciTc) has similar form to the dependence of (Afgas) on LLyafi, 
as shown in Fig. |4{e). This explains the break in this relation at 
z = 6.6, which corresponds to the switch from quiescent to burst 
galaxies, as was the case for the cold gas mass. 

Finally, the metallicity (Zcom) correlates strongly with 
Lhyafi but fairly weakly with redshift, as shown in Fig.|4{iD- The 
metallicity of galaxies with LLya [erg h~^] > lO*'^ is found to 
be around ~ 10~^. 

The quantities shown in Fig. |4] are fed into the Monte Carlo 



radiative transfer code to calculate the Lya escape fraction fcsc and 
the line profile. The net Lya luminosity of the galaxy is then simply 
Li^ya ~ fcsc^Lya.o- The value of fesc depends, of course, on the 
outflow geometry. 

4.2 Choosing the outflow parameters 

In the following, we outline the procedure to choose the value of 
the free parameters of our outflow models. Our strategy to set the 
value of these parameters ( [Jm , fv , fn] for the Shell geometry, and 
[fv , fR.] for the Wind geometry) consists of matching the observed 
cumulative Lya luminosity functions (CLFs) in the redshift range 
< z < 7. Then, we use the values obtained to make the predic- 
tions in the remainder of the paper. 

A similar strategy was followed in previous modelling of 
Lya emitters JEe Delliou et al. 2005, 2006; Nagamine et al. 200^; 
iKobavashi. Totani & Nagashimal2007i) . However, it is worth point- 
ing out that, to date, the model presented here is the only one that 
attempts to match the recently observed abundances of Lya emit- 
ters at 2; = 0.2 and at higher redshifts at the same time. 

As a first step towards our definitive model, we show in the left 
panel of Fig.|5]the cumulative luminosity function (CLF) predicted 
by GALFORM alone and applying a fixed fcsc, i.e. without using the 
Monte Carlo radiative transfer code. 

The CLF constructed using the intrinsic Lya luminosity of 
galaxies greatly overpredicts the observed estimates in the redshift 
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Figure 5. The cumulative luminosity function of Lya emitters at redshifts z = 0.2, 3.0, 5.7 and 6.6. The green lines in the left panel show the CLF obtained 
using the intrinsic (without attenuation) Lya luminosity at different redshifts. The cyan lines show the CLF obtained when applying a fixed global escape 
fraction of f^sc = 0.02 to all galaxies. The magenta dashed lines show the effect of applying an escape fraction that varies with redshift, to match the 
observational data. The values chosen are written in the labels. The right panel shows with solid red (blue) lines the CLFs obtained with the Wind (Shell) 
geometry when choosing the free parameters to match the observational CLF at z = 3, and then the same parameters were then used to predict the CLF at the 
other redshifts (see the text for details). In both panels, the observational CLFs at each redshift ai'e shown with gray symbols: At z = 0.2 d ata is taken from 
Deharveng et al. (2008) and Cowie. Barser & Hu (2 oiol): at2 = 3.0 from|Gronwall et al. (2007); Ouchi et al. (2008) and Ranch et all fcOOSh : at ^ = 5.7 from 
Shimasaku et al.l <200ei) : IOuchi et al.l <20O8D and lHu etal.U2O10l) : and at z = 6.6 from lKashikawa et al.l <2006l) and lHu et aU<2O10i) . 



range 0.2 < z < 6.6. This is not surprising, since, as discussed 
earlier, Lya photons are expected to suffer an important attenuation 
due to the presence of dust, and therefore, to have small escape 
fra ctions. Fig. [5] shows also the resul ts of the approac h followed 
by iLe Delhou et alj ( I2005L |2006|) and lOrsi et alTl2008h . in which 
fesc = 0.02 is adopted for all galaxies, regardless of their physical 
properties or redshift. This method is equivalent to a global shift 
in the intrinsic Lya CLF faintwards. This assumption is able to 
match remarkably well the observed CLF of Lya emitters atz = 3. 
Also, it is found to provide a good fit to observational estimates at 
z — 5.7 and to slightly underpredict the observed CLF at z = 6.6. 
However, the largest difference occurs at 2 — 0.2, where the fixed 
fcsc = 0.02 scenario overpredicts the observational CLF by a factor 
of ~ 5 in Lya luminosity. Note that the 2 = 0.2 CLF estimate was 
not available at the time Le Delliou et al. studied Lya emitters with 
a constant value for fosc, and therefore these authors were not aware 
of this disagreement. 

Fig. |5] also shows the effect of choosing a value for the Lya 
escape fraction that varies with redshift. At z = 0.2, we find that 
a value of fcsc = 0.005 is needed to match the observational data. 
At z = 6.6, a value of fcsc = 0.03 provides a better fit to the 
observational data than fcsc = 0.02. 

This second method, i.e. varying fcsc with redshift to find 
the best fitting value, has been also used in previous works (e.g. 
iNagamine et a l. 2010). Although it reproduces the observed abun- 
dances of Lya emitters at different redshifts, it lacks physical mo- 
tivation. Therefore, we now turn to implementing our Monte Carlo 
radiative transfer model to compute fcsc . 

In our implementation of the Lya radiative transfer code in 
GALFORM, the code can take anything from a few seconds up to 



minutes to run for a single configuration. Hence, the task of finding 
the best combination of parameters that match the observational 
CLFs shown in Fig. |5] could be computationally infeasible if we 
were to run the code on each galaxy from a GALFORM run (which 
could be in total a few hundred thousand or even millions of galax- 
ies). To tackle this practical issue we construct a grid of configura- 
tions for a particular choice of parameters values. This grid spans 
the values of the physical properties predicted by GALFORM that 
are relevant to constructing the outflows. The grid is constructed in 
such way that the number of grid points is significantly smaller (a 
factor 10 or more) than the number of galaxies used to construct it. 
An efficient way to construct the grid allows us to run the Monte 
Carlo radiative transfer code over each of the grid points, hence 
obtaining a value for fcsc in a reasonable time. Finally, we interpo- 
late the values of fcsc in this multidimensional grid for each galaxy 
in GALFORM, according to their physical parameter values, thus 
obtaining a value of fcsc for each galaxy. The methodology to con- 
struct the grid and to interpolate the value of fcsc for each galaxy is 
described in detail in Appendix Icl 

Ideally, we would like to find a single set of parameter values 
to reproduce the observed CLFs at all redshifts. Fig. |5] shows the 
result of choosing the best combination of parameters to match the 
observational CLF at 2 = 3. We chose this particular redshift be- 
cause the observational measurements of the CLF at this redshift 
span a broad range of Lya luminosities, and the scatter between 
the different observational estimates of the CLF is smaller than at 
other redshifts. 

For simplicity, we have chosen fv ~ 1.0 in both geometries. 
Also, = 0.1 was used for the Shell geometry. Hence, the only 
truly adjustable parameter in each model is /r. We have found that 
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Figure 6. The cumulative luminosity function of Lya emitters at redshifts z = 0.2, 3.0, 5.7 and 6.6. The CLF predicted using the Wind geometry is shown 
with a solid red curve, and the Shell geometry is shown in blue. The contribution of quiescent and starbursts to the model CLFs are shown with dashed 
and dot-dashed curves, respectively. Note that at z^,3 starbursts completely dominate the total CLF. Observational CLFs at each redshift are shown with gray 
symbols, like in Fig. [5] 



a value of Jr = 0.85 for the Shell geometry, and fn = 0.15 for 
the Wind geometry are needed to match the observed CLF of Lya 
emitters at 2: = 3. 

Despite finding a reasonable fit to the observational data at 
this redshift, neither of our outflow geometries is able to match the 
observed CLFs over the whole redshift range < 2: < 7 if one 
fixed set of parameters is used. 

As shown in Fig.O the observed CLF at z = 0.2 is particu- 
larly difficult to match when the model parameters have been set to 
match the CLFs at higher redshift. This occurs because the num- 
ber density of Lyg emitters at z = 0.2 reported observationally 
jPeharveng et alJl2008l ; ICowie. Barger & Hij|2010h is much lower 
than what our model predicts, implying feac is also considerably 
lower than what our model suggests. On the other hand, the ob- 
served CLF at 2: = 6.6 implies a higher abundance of Lya emit- 
ters than what our model for computing the Lya escape fraction 
predicts. 

An improved strategy to reproduce the observed CLFs of Lya 
emitters at < z < 7 is to choose the free parameters of the 
models in the following way: 

(i) The expansion velocity of the outflows is set to be equal to 
the Lya-weighted circular velocity of the galaxies (i.e. fv~lm 
both geometries). 



(ii) In the Shell geometry, the fraction of cold gas mass in the 
outflow is fixed at /m = 0.1. 

(iii) Given the strong evolution of the half-mass radius of 
galaxies with redshift, the value of /j? is allowed to evolve with 
redshift. 

Galaxies at low redshift are predicted to be larger in size than 
galaxies at higher redshifts, so if fn is fixed like the other parame- 
ters, then outflows at lower redshift would have, on average, larger 
sizes than at high redshifts, making the associated column densities 
smaller (see Eqs.[T]and|19t. If the other outflow properties do not 
evolve strongly with redshift, then galaxies at low redshift would 
have larger Lya escape fractions than at higher redshifts. As dis- 
cussed above, this is the opposite trend needed to reproduce the 
observed CLFs. 

However, there is growing evidence that star forming re- 
gions in local ultraluminous infrared galaxies (ULIRGS) are sig- 
nificantly smaller than similarly luminous ULI RGS and submil- 
limetre galaxies (SMG s) at higher redshifts (see llono et al]|2009l : 
iRuiopakarn et alj|201 1 for a comparison of sizes). If the outflow 
radius in starbursts is assumed to scale with the size of the star 
forming regions instead of the full galaxy size, then fn would have 
a natural redshift dependence in the direction we need. 

Hence, we employ a simple phenomenological evolution of 
fn with redshift. We call f^i the radius parameter of starbursts and 
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Table 1. Summary of the parameter values of the Shell and Wind geometries 
used to fit the Lya cumulative luminosity function at different redshifts (see 
the text for details). 



allow it to scale with redshift like a power-law: 



(21) 



where g and 7 are free parameters. Since there is no equivalent 
observational evidence for the size of star forming regions in quies- 
cent galaxies scaling with redshift, we set f^, the radius parameter 
for quiescent galaxies, to be an adjustable parameter but fixed (i.e. 
independent of redshift). 

Table[T]summarizes a suitable choice of the parameters values 
used in our model. Also, Fig. [6] shows the predicted cumulative 
luminosity functions obtained with that choice of parameters. 

Our previous modelling of Lya emitters used the simple as- 
sumption of a constant Lya escape fraction, with fosc ~ 0.02 be- 
ing a suitable v alue to reprod uce the Lya CL Fs at 3 < z < 7 
( iLe Delliou e~ al. 2005, 2006; Orsi et alj[20o3) . As shown in Fig. 
[S] this simple model overestimates the CLF at z = 0.2, but it re- 
produces remarkably well the CLFs at higher redshifts. Our out- 
flow geometries, on the other hand, are consistent with the obser- 
vational CLFs at all redshifts, although they fail to reproduce their 
full shape. This may be surprising at first, since the intrinsic Lya 
CLF (shown in green in Fig.|5]( roughly reproduces the shape of the 
observed CLFs, although displaced to brighter luminosities (which 
is why the constant escape fraction scenario works well at repro- 
ducing the CLFs). However, in our model, the escape fraction in 
each galaxy is the result of a complex interplay between several 
physical properties, and this in turn modifies the resulting shape of 
the CLF. 

The sizes of the outflows predicted by our model, com- 
pared to the extent of the galaxies themselves (quantified by the 
half-mass radius {R1/2)) are very different between the two out- 
flow geometries. In quiescent galaxies these are 1.5 and 20 per- 
cent of the half-mass radius of the galaxies, in the Wind and 
Shell geometries respectively. Similarly, at 2; = 0.2, outflows 
in starbursts are 2 and 26 percent of the half-mass radius in 
the Wind and Shell geometries. The rather small size of out- 
flows in the Wind geometry appears to be in contradiction with 
observations of Lya in local starbursts which display galactic- 
scale outflows (see, e .g. iGiav alisco. Koratkar & Calzetti 1996; 
Thuan&Izotov 1997; Kunth et alj|1998l; IMas-Hesse et alj 120031: 



OstUn et al. 2009 : Mas-Hesse et alj2009h . However, local starburst 



samples are sparse and still probably not large enough to charac- 
terise the nature (in a statistical sense) of Lya emitters at low red- 
shifts. 

At higher redshifts, our model keeps the sizes of outflows in 
quiescent galaxies unchanged (with respect to their half-mass ra- 
dius). However, outflows in starbursts grow in radius, relative to 
their host galaxy, according to a power law, as given by Eq. l l21t . 
with the best fitting values listed in Table [T] so at 2; = 3 their sizes 
are 27 and 80 percent of the half-mass radius for the Wind and 
Shell geometries, respectively. By z = 6.6, the sizes are 110 and 
145 percent of the half-mass radius. Therefore, at 2; > 3, all out- 
flows in Lya-emitting starbursts are galactic-scale according to our 
models. 
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Figure 7. The neutral hydrogen column density of outflows for galaxies as 
a function of Lya luminosity at redshifts z = 0.2 (top), z = 3.0 (middle) 
and z = 6.6 (bottom). The Wind geometry is shown in red (left), and the 
Shell geometry in blue (right). Filled circles show the median of the column 
density distribution as a function of the attenuated Lyo luminosity. Open 
circles show the same but as a function of the intrinsic Lyo luminosity. 
Error bars show the 10-90 percentiles of the column density distribution at 
each Lya luminosity bin. 



A consequence of setting the free parameters in the model to 
reproduce the low-z data is that the CLF of Lya emitters at z — 0.2 
has a contribution from both quiescent and starburst galaxies. De- 
spite the details over which component is dominant at this redshift 
(which is something somewhat arbitrary given the freedom to ad- 
just the other free parameters of the models), the CLF of Lya emit- 
ters at high redshifts is invariably dominated by starbursts, and only 
a negligible fraction of Lya emitters are quiescent galaxies. 



5 PROPERTIES OF Lya EMISSION 

Having chosen the parameters in our outflow geometries we pro- 
ceed to study the predictions of our hybrid model for the properties 
of Lya emitters. Whenever possible, we compare our predictions 
with available observational data. 



5.1 Column densities 

An immediate consequence of the choice of parameters shown in 
Table[T]is the distribution of predicted column densities of the out- 
flows of galaxies. The derived hydrogen column density is shown 
in Fig. [T] as a function of Lya luminosity for different redshifts. 
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Figure 8. The escape fraction as a function of Lya luminosity at ^ = 0.2 
(top), 2 = 3.0 (middle) and z = 6.6 (bottom). The Wind geometry predic- 
tions are shown in red (left), and the Shell geometry in blue (right). Solid 
circles show the median of the escape fraction at different attenuated Lya 
luminosity bins. Open circles show the same relation as a function of in- 
trinsic Lyo luminosity instead. Error bars represent the 10-90 percentiles 
around the median of the escape fraction distribution at a given Lya lumi- 
nosity bin. 

Here we study this distribution both as a function of intrinsic and 
attenuated Lya luminosity. 

When studied as a function of intrinsic Lya luminosity, our 
model displays a large range of column densities, varying from 
Nh ~ 10^^ [cm^^] to Nh ~ 10^'* [cm^^], depending on the out- 
flow geometry and the redshift. Typically, in both outflow geome- 
tries, the column density increases with intrinsic Lya luminosity, 
which is due to the overall increase in the mass of cold gas and the 
mass ejection rate with intrinsic Lya luminosity. However, in some 
cases, there is a noticeable decrease of the column density with in- 
creasing intrinsic Lya luminosity, reflecting the rather complicated 
relation between the quantities which affect the column density and 
the intrinsic Lya luminosity, as shown in Fig.|4] 

When including attenuation due to dust, this relation is mod- 
ified, as shown by the solid circles in Fig. |7] The number of scat- 
terings scales with the column density of the medium, hence the 
escape fraction is low (or effectively zero in some cases, accord- 
ing to our model) for outflows with large column densities. As a 
result, the column density distribution as a function of attenuated 
Lya luminosity spans Nh ~ 10^^""^^ [cm~^]. Also, galaxies with 
brighter (attenuated) Lya luminosities tend to have smaller column 
densities. 

The column densities predicted by our models are similar to 
those inferred bv lVerhamme et al] j2008l) on fitting their models to 



a small sample of high redshift galaxies with high resolution spec- 
tra. These authors fit Lya line profiles with a Monte Carlo radia- 
tive transfer model using a geometry identical to our Shell geom- 
etry. In addition, Verhamme et al. present a compilation of results 
for observationally-measured column densities of local starbursts 
showing Lya emission, with values between 10^** — 10^^[cm~^], 
which are also consistent with our model predictions. 



5.2 Lya escape fractions 

A fundamental prediction of our models is the distribution of Lya 
escape fractions, shown in Fig. [8] In terms of the attenuated Lya 
emission, we find that fosc grows monotonically from 10""^ in 
the faintest Lya emitters (with Z/Lya ~ 10^*[ergs-ih-2]) to 
fcsc ~ 0.3 at I/Lya ~ 10'*" [erg ^ h^^] . Brighter galaxies have 
in general escape fractions ranging from 0.1 to ~ 1, depending on 
the redshift and outflow geometry. Perhaps not surprisingly, both 
models predict similar distributions of escape fractions as a func- 
tion of Lya luminosity, since the models are forced to match the 
observed CLFs (Section r4.2t . 

On the other hand, in terms of the intrinsic Lya emission. Fig. 
[8]shows that only in some cases it is possible to obtain a median fcsc 
above zero. This reflects the complicated interplay of physical con- 
ditions which shape the escape fraction of Lya photons. In other 
words, the value of the intrinsic Lya luminosity in a galaxy does 
not determine fcsc and, therefore, its attenuated Lya luminosity. 

Fig. [8] illustrates further the contrast betweetn using a phys- 
ical model to compute the Lya escape fraction, and the constant 
fcsc scenario. In the former, galaxies seen with high Lya luminosi- 
ties have high escape fractions, so their intrinsic Lya luminosities 
are similar to the attenuated ones. Fainter galaxies, on the other 
hand, have lower escape fractions, meaning that their intrinsic Lya 
luminosity is orders of magnitude higher than their observed lumi- 
nosity. This implies that a Lya emitter with a given intrinsic Lya 
luminosity could have either a high or low value of fcsc, depending 
on its physical characteristics (galaxy size, metallicity, circular ve- 
locity, SFR, etc.). In the constant escape fraction scenario, on the 
other hand, all galaxies at a given observed Lya luminosity have 
the same intrinsic luminosity, which means that galaxies need only 
to have a high intrinsic Lya luminosity to be observed as Lya 
emitters. 

We conclude from Fig.[8]that a high production of Lya pho- 
tons (or, equivalently, high intrinsic Lyman continuum luminosity) 
does not guarantee that a galaxy is visible in Lya. Our modelling 
of the escape fraction selects a particular population of galaxies to 
be observed as Lya emitters, which in general is found to have 
lower metallicities, lower SFRs and larger half-mass radii than the 
bulk of the galaxy population. This is discussed in more detail in 
Section|6] 

It is worth asking at this point whether the predicted escape 
fractions are consistent with observational estimates. Observation- 
ally, fcsc is generally calculated either by inferring the SFRs from 
the Lya luminosity and comparing to the SFR estimated from 
the UV coritinuum (e.g. iGawiser et al] l2006l : Isianc et al.l |20I ll: 
iGuaita et al.ll20Io l). or by using the ratio between Lya and another 
non-resonant hydrogen recombination line. The former requires as- 
sumptions about the stellar evolution model, the choice of the IMF 
and the modelling of dust extinction. The second method only re- 
lies on the assumption that the comparison line is not affected by 
resonant scattering, and thus the extinction can be estimated reli- 
ably from an extinction curve, and that the intrinsic line ratio cor- 
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Figure 9. The Lyo escape fraction as a function of tlie observed (i.e. at- 
tenuated) ratio of Ha to H/3 flux at 2 = 0.2. Red circles sliow tlie median 
predictions using tlie Wind geometry, and blue circles show the median 
predictions of the Shell geometry. En'or bars in both cases show the 10-90 
percentil es of the distributi on. Grey symbols show observational measure- 
ments bv lAtek et ai] j2008h (crosses), and GALEX and lUE samples (dia- 
monds and triangles respectively) from lAtek et alj i2009l) . 



responds to the ratio of the emission coefficients, assuming case B 
recombination. 

For this reason, we focus on the escape fraction measured us- 
ing line ratios. These are often presented as a function of colour 
excess E{B — V), which in turn is estimated from non-resonant 
recombination lines such as Ha and H/3 together with an assumed 
extinction curve, which describes foreground extinct ion jAtek et al] 
I2OO8L l2009l : lOstUn et alJ l2009l : iHaves et alJ I2OI0I) . For the pur- 
poses of comparing the model predictions to the observationally- 
estimated values of fesc , we convert the values of E{B~ V) quoted 
in the observations to Rohs = /Ha,obs//H/3,oba, the ratio between 
the observed fluxes of Ha and H/3, since this ratio is a direct pre- 
diction from G ALFORM. To compute i?obs, we follow a standard 
relation jAtek et al. 2008 ). 



l0g(i?obs) = l0g(i?int) - E{B - V) 



fc(A,)^fc(A^) 
2.5 



(22) 



where Rint = 2.86 is the intrinsic line ratio between Ha and H/3 
typically assumed under Case B recombination fo r a medium at 
a temperature of T = 10000[K] jOsterbrock|[T989l) . and k{Xa) = 
2.63, and k{Xp) = 3.71 are the values of the normalized extinction 
cu rve at each corresponding wavelen gth from the extinction curve 
of lCardelli. Clavton & Ma"ai3 i ll989l) . 

Fig. |9] shows the predicted relation between the Lya es- 
cape fraction and Rohs compared to observational estimates from 
lAteketal.U200 ^ and an analys is of UV spectros copic data from 
the GALEX and lUE surveys bv lAtek et alTj2009l) . 

The model predictions shown in Fig. [9] include only galaxies 
with log(LLj;a [erg s^^ h^^]) > 41.5, in order to approximately 
reproduce the selection of Lya emitters in the GALEX sample. 
Note, however, that the observational points shown in Fig.|9]do not 
represent a complete statistical sample. Therefore it is not possible 



to perform a fair comparison, and the results shown here should be 
regarded as illustrative. 

The Shell geometry shows remarkable agreement with the ob- 
servational estimates of fcsc , reproducing the trend of lower escape 
fractions in galaxies with larger i?obs- The Wind geometry, on the 
other hand, is only partially consistent with the observational data, 
and it predicts a rather flat relation between fesc and i?obs- 

5.3 The Lya Equivalent width distribution 

The equivalent width EW measures the strength of the line with 
respect to the continuum around it. We compute the EW simply by 
taking the ratio of the predicted Lya luminosity of galaxies and the 
stellar continuum around the Lya line as computed by GALFORM, 
including attenuation by dust. Fig. [TO] shows a comparison of the 
EW distribution measured at different redshifts with the predictions 
from our outflow geometries. 

Overall, both outflow geometries predict EW distributions 
broader than the observational samples for all redshifts studied. To 
characterise the predicted EWs, we compute the median of the dis- 
tributions. This will depend on the Lya luminosity limit applied 
to the sample to make a fair comparison with observed data. In 
Fig. llO | we compare ou r model predictions w ith obser vational data 



from 



Cowie. Barger & Hu (2010) at 2 
5.7. 



: 0.2. IOuchie tal. (2008)at 



z — 3.1 and also at z 

The comparison between the observed and predicted EW by 
the Wind geometry median values (shown by vertical dashed lines 
in the left panel of Fig. 1101) is encouraging at 2 = 0.2, where the 
EW distributions have a median value of ~ 30 A. However, both 
geometries give consistently higher median values of the EW when 
compared to the observational data at higher redshift. In addition, 
the disagreement in the median values of the EW distribution be- 
tween observational data and the Shell geometry predictions be- 
comes larger as we go to higher redshifts. 

At 2 = 0.2, the EW distribution predicted using both 
outflow geometries is broa der than the observational sample of 
ICowie, Barger & Hul ilOld) and reaches values of EW ^ 300A, 
whereas the observational sample only reaches EW ~ 130A. 

At 2 = 3.0, the observational sample seems to peak around an 
EW « 8OA and then declines until reaching EW ^ 150 A. The 
predicted EWs are consistent with the observational distributions, 
although the former reach values as high as EW ~ 400 A. A 
similar disagreement between model predictions and observational 
data is found at 2 = 5.7 

The disagreement found between the EW values in the model 
predictions and the observations is difficult to understand from 
studying only the EW distributions. Therefore, we perform a more 
detailed comparison by studying the relation between the median 
EW and the Lya luminosity, as shown in the right panels of Fig. 

m 

Fig.[TO]shows that there is good agreement between the model 
predictions and the observational data at z — 0.2, where both 
models seem to reproduce the range and scatter of the observa- 
tions. Overall, both models predict an increase of EW towards 
brighter Lya emiters: The Wind geometry predicts that the me- 
dian EW increases from « 30A to ~ lOOA for Lya luminosities 
Lhya ~ 10*" — 10*^ [erg whereas the Shell geometry 

predicts a much steeper increase of the median EWs, from ~ 3A to 
^ lOOA over the same Lya luminosity range. 

At 2 = 3, the EW values from the observed data are also con- 
sistent with both model predictions. The Shell geometry predicts 
a steep increase in the EW values with increasing Lya luminosity. 
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Figure 10. The (rest-frame) equivalent width distribtion at z = 0.2,3.0 and z = 5.7. The left panels show histograms of the distribution of EWs at different 
redshifts with the applied lower luminosity limit given in the boxes. The right panels show the median of the distribution of EWs as a function of Lya 
luminosity. In all cases, the Wind geometry predictions are shown in the left column in red, and the Shell geometry in the rig ht column in blue. The error 
bars on the right panel denote the 10 -90 percentiles of the distribution of EWs. Observational data (shown in gray) is taken from lCowie. Barger & h3 feOlOl) 
at z = 0.2 and lOuchi et al.l j2008l) for redshifts z = 3 and z = 5.7. The vertical dashed lines in the left panel correspond to the median values of the 
distributions. 



whereas the Wind geometry shows a rather flat relation. At z — 5.7 
the models are also consistent with most of the observed EW values 
within the predicted range of the EW distributions. 

5.4 UV continuum properties of Lya emitters 

The variant of GALFORM used in this work has been previously 
shown to match the abundance of LBGs (cha racterised by their 
UV luminosities) over a wide range of re dshifts teaugh et alj2005l : 
iLacev et alj|201 ll : iGonzalez et alj|20l" 11). Therefore, a natural pre- 
diction to study with our model is the UV LF of a Lya-selected 
sample. 

Fig. E] shows the UV (1500 A) LF of Lya emitters at red- 
shifts 3.0, 5.7 and 6.6. To compare the model predictions with ob- 
servational data, we mimic the UV selection applied to each sam- 
ple. These correspond to constraints on the limiting Lya luminosity 
and the minimum EW. Moreover, we have chosen to compare our 
model predictions with two observational samples at each redshift. 
Although the observational samples show similar limiting Lya lu- 
minosity at each redshift, the value of the minimum EW, which is 
different among observational samples, has an important impact on 
the predicted UV LF of Lya emitters, as shown in Fig.llll 

At 2 = 3, we compare our model predictions with the UV 



LF from lOuchi et alj ( l2008t) and lOronwall et all ( |2007|) . Overall, 
our models are found to underpredict the observational estimates, 
alt hough they are consist ent with the faint end of the LF measured 
bv lGronwalletal.l ( l2007t) . and in reasonable agreement at the bright 
end with the observational sample ofjOuchi et al. (2008). Likewise, 
at z = 5.7 both outflow geometries are found to undepredic t the 
UV LF of Lya em itters measured bv lShimasaku et al] ( l2006t) and 
lOuchietalJfeOOSi) . 

The situation is somewhat different at z = 6.6. Both out- 
flow geometries are con sistent with the faint end of the UV L F 
of Lya emitters shown in lKobavashi. Totani & Nagashimal ( |201(]|) . 
who recalculated the U V LF of Lya emitters from the sample of 
iKashikawa et alj (l2006l) after applying a brighter Lya luminosity 
cut, to ensure that a more complete sample was analysed. Also, 
both outflow geomet ries are consistent with the U V LF of Lya 
emitters measured bylc owie. Hu&Songaii3 l l201ll) . 

For comparison, we show in Fig.[TT]the UV LF of all galaxies 
predicted by GALFORM, i.e. without applying any selection. Note 
this was alr eady shown to agree remarkably well with observed LFs 
of LBGs in iLacevetalJ ( 1201 ll) . This total UV LF is always above 
the UV LF of Lya emitters predicted by the outflow geometries. 

At z = 3, we notice that the UV LF of all galaxies is above 
the observed LFs of Lya emitters. This is consistent with the idea 
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Figure 11. The UV luminosity function of Lyo emitters at redshifts z = 
3.0 (top), 2 = 5.7 (middle) and z = 6.6 (bottom). The Wind geometry pre- 
dictions are shown in red, and those of the Shell geo metry in blue. Obse rva- 
tional data, shown wit h gray sym bols, is taken from |Ouchi et al.l )2008h and 
iGronwall et"ai] J2007l) at 2 = S.IShimasaku et alj I2 OO6) and O uchi et alj 
i2008l ) at z = 5.7 and lKobavashi. Totani & Nagashima t2010) at z = 6.6. 
The limiting Lyo luminosity and EW used to construct the models UV LF 
are shown in each panel, in units of log{[erg h~^]) and [A], respec- 
tively. The dashed green curves show the UV LF of all galaxies, as predicted 
by GALFORM. without imposing any Lyo selection. 



that Lya emitters constitute a sub-sample of the galaxy popula- 
tion. However, at 2 = 5.7 and z = 6.6 this UV LF matches 
remarkably well the observed UV LFs of Lya emitters. A sim- 
ilar finding has also been reported in o bservational papers (e.g. 
iKashikawa et alj2006l : IOuchi et al.ll2008h when comparing the UV 
LF of LBGs to th at of Lyg emitters. Interesting ly, only the obser- 
vational LF from lCowie. Hu & Songaill JloT l|) falls significantly 
below the predicted total UV LF of all galaxies, and is at the same 
time consistent with the predicted UV LF of Lya emitters for both 
outflow geometries. 

It has been argued that attenuation by the IGM might play 
a significant role in sh aping the LF of Lya emitters at 2: > 5. 
IKashikawa etall ( I2OO6I) interpret their measured LF of Lya emit- 
ters as evidence of an abrupt decrease in the amplitude of the bright 
end of the Lya LF at 2 = 6.6 compared to 2 = 5.7. Since the 
corresponding UV LFs of Lya samples do not seem to evolve in 
the same way, this has been suggested as a result of a change in 
the ionization state of the IGM (for which Lya photons are more 
sensitive than continuum photons). 

Our model does not compute any attenuation by the IGM. The 
effect of reionization in GALFORM is modeled simply by preventing 
the cooling of gas in haloes of a given circular velocity (Kirc = 
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Figure 12. The fraction of Lya emitters with rest-frame equivalent width 
EW^i > 25A for galaxies with -19.4 < Muv - 5 \og{h) < -17.9 (as- 
terisks) and —20.9 < Muv — 51og(/i) < —19.4 (diamonds) as afunction 
of redshift. Predictions of the Wind and Shell geometries a re shown in red 
and blue, re spectively. The observ ational measurements of Sche nker et alj 
( 2011) and iPentericci et alj j201lh are shown in gray and green, respec- 
tively. 



30km/s in our GALFORM variant ), when t heir redshift is smaller 
than a redshift 2rcion ~ W (see ISenson et al. 2002, for a more 
detailed model). 

Some theoretical work exploring the abundances of Lya 
emitters at these high redshifts has attempted to take into account 
the attenuation of the Lya line due to absorption by the IGM 
(iKoba yashi, Totani & Nagashima 2007; Dayal, Mase lli & Ferraral 
l20rii ; iKobavashi, Totani & Nagashim al 120101) . In our model, a 
simple mechanism to improve the agreement between the observed 
and predicted UV LFs of Lya emitters at high redshift would be 
to add a constant attenuation factor of the Lya luminosity by the 
IGM at high redshifts, set to match the UV LF of Lya emitters. 
We defer doing this to a future paper 

We now focus on the fraction of galaxies exhibiting Lya in 
emission predicted by our models. As described in Appendix [Cl in 
our radiative transfer model we follow a maximum of 1000 pho- 
tons per galaxy, and so we can compute a minimum value for the 
escape fraction of 10~^ . However, in a significant fraction of galax- 
ies, none of the photons escape from the outflows, resulting in them 
being assigned fcsc ~ 0, and thus having no Lya emission at all. 
These galaxies could be related to th e observed populatio n of galax- 
ies showing Lya in absorption (e.g. IShaplev et alj2003h , which we 
will examine in a future paper. 

Fig. [T2] compares our model predictions with the observed 
fraction of Lya e mitters found in Lyma n-break galaxies (LBGs) 
at 2 ^ 4 - 7 by ISchenker eTal] j201lh and by IPentericci et alj 

ij) at 2 ~ 7. Lya emitters are defined here as galaxies with a 
Lya EWrf ^ 25A. The samples are split according to two dif- 
ferent rest-frame UV magnitudes ranges, as shown in Fig. [12] For 
simplicity, we define UV magnitudes at rest-frame 1500A. 

For the two ranges of UV-magnitudes shown, the Shell geom- 
etry appears to have a larger fraction of Lya emitters as a function 
of redshift than the Wind geometry. W hen comparing to the ob- 
servational data of iSchenker et alj ( l201lh . we find that overall both 
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Figure 13. Composite Lyo line profiles of s amples at different redshi fts as indicated in eac h panel. The gray shaded regions show the Lya line profiles from 
composite observational samples taken from lOuchi et al.l bOOSlEoioh and lHu et al.l <201(ll) (as indicated in the legend of each box), whereas the solid curves 
show the model predictions using the Wind geometry (red) and the Shell geometry (blue). Galaxies in the models used to construct the composite spectra have 
been selected in the same way as the observed spectra, with the criteria shown in the legend of each box, where Li^ya is shown in units of [erg h~^] and 
the EW in A. 



outflow geometries underpredict the fraction of Lya emitters mea- 
sured observationafly. In particular, only at 2; > 5 are the Wind and 
Shell geometries consistent with the fraction of Ly« emitters mea- 
sured in the UV magnitude range —20.9 < Mjjv ^ 51og(/i) < 
-19.4. 

Despite the differences, both outflow geometries predict an 
increase in the fraction of Lya emitters with redshift, which is 
qualitatively consistent with the observati ons up to z ^ 6. A t 
hig her redshifts, the obser vational data of ISchenker et alj ( l201lh 
and IPentericci et alj ( I2OIII) suggests a decline in the fraction of 
Lya emitters at z ~ 7. They interpret this decline as the impact 
of the neutral IGM attenuating the Lya luminosity from galaxies. 

Our models, on the other hand, show a trend consistent to the 
observations in the magnitude range —19.4 < Muv — 5 \og{h) < 
— 17.9, except at the highest redshift, z = 7.3 where both models 
increase their fraction of Lya emitters instead of decreasing it, as 
observations do. However, it is worth noticing that in the UV mag- 
nitude range —20.9 < Muv — 51og(ft) < —19.4 both outflow 
geometries predict a decline of the Lya fraction similar to the one 
found observationally. Therefore, Fig.ll2lshows that our model pre- 
dictions imply that the decline in the Lya fraction in LBG samples 
found at high redshifts is not conclusively driven by the presence 
of neutral HI in the IGM attenuating Lya photons. 



5.5 Observed Lya line profiles 

Observational measurements of individual and stacked line pro- 
files of Lya sugge s t the presence of outflows in galaxies (e.g . 
Shaplev et al] l2003l: iK ashikawa et al.' '200^: 'Pawson et al .l |2007|: 
" IHu et al. 2010; Kornei et al. 2010; Ouc hi et alj 
12010 ; IKulas et al.ll2011i') . Lya emitters can be 



Ouchi et alj|2008l; 
20ld ;'Steidel etal 



characterised by studying the spectral features of the composite 
spectrum from a set of spectroscopic observations. The most promi- 
nent feature observed are asymmetric peaks, where the line is ex- 
tended towards the red side. Other common spectral features are 
the appearanc e of a secondary pe ak and P-Cygni absorption fea- 
tures (see, e.g IShaplev et al]|2003l) . 

In this section we compare our model predictions with the 
comp osite spectra of high redshift samples of L ya emitters stud - 
ied bv lOuchi et id] I2OO8), Ouchi et al. ( 2010) and lHu et alj ilOld) . 
Based on the method used to construct composite spectra in the ob- 
servational studies, we construct composite spectra in the model as 
follows. First, the Lya profiles are normalised to their peak values. 
Then the spectrum is shifted so that the peaks coincide with the 
Lya line centre, Alyq = 1216A. Finally, the spectra are averaged 
at each wavelength bin. 

The above method for constructing a composite spectrum has 
some important drawbacks. Since the redshift of Lya emitters is 
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computed from the wavelength of the peak of the Lya hne, any 
offset of the peak due to radiative transfer effects is removed (see 
Fig. [Hand Fig. [2}. On the other hand, the normalisation of the line 
profiles to the peak value can help enhance certain spectral features 
inherent to Lya emitters by removing any dependence of the com- 
posite spectrum on Lya luminosity. This also means that spectral 
features characteristic of a particular Lya luminosity could be dif- 
ficult to spot. 

The top panels in Fig. [T3] show a comparison be tween com- 

[ osite spectr a at redshifts z = 3, 5.7 and z = 6.6 from lOuchi et alj 
2OO8 II2OIOI) with the predictions from both outflow geometries. 
Overall, both outflow geometries show very similar composite line 
profiles, regardless of redshift or limiting luminosity. A hint of a 
secondary peak redward of the line centre is weakly displayed for 
some configurations, but it is not strong enough to make a clear 
distinction between the line profiles predicted by both geometries. 

The Lya line in the composite spectrum estimated at 2: = 3.1 
bv lOuchi et al.l ( l2008i) is broader than our model predictions. This 
could suggest that these galaxies have a larger column density or 
expansion velocity than our model predicts, as Fig. |2] also shows. 
A\. z = 5.7 and z — 6.6, Fig. 1131 shows remarkable agreem ent 
between the models and the observations bv lOuchietalldioOSl) and 
lOuchietaljJioToh . respectively. 

The situation is different when compar ing with the com posite 
spectra at z = 5.7 and 2 — 6.6 estimated bv lHu et al] ( l2010 l). Their 
composite spectrum is constructed from samples of Lya emitters 
of similar limiting Lya luminosity but with EWrt > lOOA , signif- 
icantly greater than the EW limit in the lOuchiet"ai] ( l2008h samples 
at 2 > 5, which have (EWrt > 20A). 

In this case, the observational composite spectra are narrower 
than our model predictions. Moreover, the observed sp ectra appear 
to be more asymmetric than their counterparts in the lOuchi et alj 
( l2008h and lOuchi et al. IIZOl^) samples. The asymmetry of the Lya 
line varies with the outflow expansion velocity, as discussed in Fig. 

H 



It is interesting to notice that the composite spectra o f lHu et alj 

( l2010h show a sharp cut-off on their blue side, which is not well 
reproduced by our outflow geometries. This lack of photons in 
the blue side of the spectra could be interpreted as the impact 
of the IGM at these high red shifts removing the blue-side o f the 
Lya spectrum, as discussed bvlDiikstra. Lidz & Wvithd (120071) an d 
more recently by iLaursen. Sommer-Larsen & Razouniiov I2OIII) . 
Since our outflow geometries do not show this feature, it seems 
likely to be caused by the presence of a neutral IGM. However, this 
does n ot explain the overall good agreement with the Ouchi et al] 
( |2008 |) composite spectra at 2 = 5.7 and the lOuchi et al. 1 2010l) 
composite spectra at 2 = 6.6. 

The predicted Lya profiles shown in Fig. [13] show reasonable 
agreement with the data, thus supporting the idea that Lya photons 
escape mainly through galactic outflows. It is worth reminding the 
reader that our outflow geometries are not tuned to reproduce the 
observed line shapes, and therefore these represent genuine predic- 
tions of the outflow geometries. 



6 THE NATURE OF Lya EMITTERS 

After performing the detailed comparison between observational 
data and model predictions in the previous section, we conclude 
that our outflow geometries reproduce at some level the general 
features of the observations. We now turn to the question of what 
are the physical conditions that make a galaxy observable through 



its Lya emission. As previously shown in Fig. [12] only a fraction 
of (JV^-selected galaxies have detectable Lya emission, implying 
that this selection targets galaxies with particular characteristics. To 
reveal the properties of Lya-selected galaxies, we show in Fig. 1 141 
a comparison between all galaxies and what we define here as a 
typical Lya emitter, i.e. a galaxy with log{LLya [erg h~^]) ^ 
41.5 and EWii ^ 20A. This corresponds to a Lya luminosity limit 
where Lya emitters are abundant over the redshift range 0.2 < 
z < 6.6, and the EW limit corresponds to a typical EW limit in 
observational samples. We use a top-hat filter centered on a rest- 
frame wavelength of A = 1500A to ensure that the rest frame UV 
magnitude is the same for all redshifts. 

The top row of Fig. [14] shows that, in general, both outflow 
geometries predict Lya emitters at high redshift to have similar or 
somewhat lower metallicities than the bulk of the galaxy population 
at the same UV magnitude. At 2 = 0.2, however, the metallicity 
of Lya emitters is larger than that of the overall galaxy population 
for a range of UV magnitudes. 

The result that Lya emitters have lower metallicities at higher 
redshifts than the overall galaxy population may not be surprising, 
since to first order we expect galaxies with low metallicities to have 
a low amount of dust and thus to be less attenuated in Lya emis- 
sion. However, according to Eq. illh metallicity is not the only 
factor controlling the amount of dust, which is why our models 
predict that at some magnitudes, typical Lya emitters can have the 
same (or even higher) metallicities than the bulk of the galaxy pop- 
ulation. 

Observational studies of the Lya emitter population at 2 ~ 
0.2 have found that these galaxies have in general lower metal- 
licities than the bulk of the galaxy populati on at the same 
stella r masses instead of UV m agnitude (e.g. iFinkelstein et alj 
l201lt ICowie. Barger & Hul l2oToh . At higher redshifts the same 
conclusion is drawn from the observational samples (see, e.g. 
Gawiser et al. 2006; Pentericci et al. 2007; Finkelstein et al. 200^. 
The metallicities our outflow geometries predict at different stellar 
masses cannot be directly compared with observational estimates, 
since these compute the stellar mass from SED fitting assuming a 
Salpeter or similar IMF. Our model, on the other hand, assumes a 
top-heavy IMF in starbursts, which are the dominant component of 
the Lya emitter population in our model (see Fig.|6j. 

Historically, it was thought that metallicity was the main 
parameter driving the observability of the Lya line, imply- 
ing that Lya emitters shou l d be essentially metal-free galax- 
ies jMeier & Terl evichI | l98l| ; iHartmann. Huchra & Gelle3 119841; 
iHartmann et alj 1988). Our model predictions, on the other hand, 
imply that the observability of the Lya line depends on the inter- 
play between more physical properties. Nevertheless, it is not the 
metallicity itself that determines the amount of dust that photons 
need to cross through when escaping. In our models, the dust op- 
tical depth depends both on the metallicity of the galaxies and on 
the hydrogen column density of the outflows, which is in turn a 
function of several properties, such as size, mass ejection rate, ex- 
pansion velocity and cold gas content. Therefore, it is not surprising 
that Lya emitters can be either metal-poor or metal-rich compared 
to the bulk of the galaxy population. 

Fig. [14] also reveals that typical Lya emitters in the Wind ge- 
ometry have overall lower metallicities than their counterparts in 
the Shell geometry. This result is consistent with the higher sensi- 
tivity to dust in the Shell geometry compared with the Wind geom- 
etry, shown in Fig.[T] Due to the different response to dust of each 
model, it is natural that the predicted metallicities are correspond- 
ingly different. 
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Figure 14. Galaxy properties as a function of extincted UV magnitude for redshifts z = 0.2 (left), 2 = 3.0 (middle) and z = 6.6 (right). Each row shows 
a different property. The top row shows the weighted metalUcity of the cold gas, the middle row the instantaneous SFR, and the bottom row the half-mass 
radius. Green circles denote the median of the distribution including all galaxies per magnitude bin. Red and blue circles denote the median of the distribution 
of galaxies selected as typical Lya emitters (i.e. with log(LLj^Q [erg h~^]) ^ 41.5 and EW 20 A) in the Wind and Shell geometries, respectively. 
Error bars denote the 10-90 percentiles of the corresponding distribution. 



Another clear difference between the overall galaxy popula- 
tion and Lya emitters is seen when looking at the instantaneous 
star formation rate at different UV magnitudes, shown in the mid- 
dle panel of Fig. [14] Here it is clear that regardless of redshift or 
UV magnitude, both models predict that Lya emitters should have 
smaller SFRs than typical galaxies with the same UV magnitude. 

The instantaneous SFR has an important impact on the escape 
of Lya photons in both models. In the Wind geometry, the mass 
ejection rate of the outflow depends directly on the SFR, as is clear 
from Eq. J13l >. The mass ejection rate is in turn directly propor- 
tional to the hydrogen column density of the outflows, according to 
Eq. ( 1191 ), which in turn affects the path length of photons, the num- 
ber of scatterings, and the amount of dust, as discussed previously. 



Hence, a galaxy with a low SFR could have a high Lya escape 
fraction, making it observable. In the Shell geometry the SFR is 
also important, since it correlates with the cold gas mass, which in 
turn determines the hydrogen column density of the outflow. 

Our outflow geometries are consistent with the observational 
evid ence for Lya emit t ers having modest star formation rate s 
(e.g. lGawiser et all 20061 : lOronwall et al.ll200llGuaita et alj201(]h . 
However, it is worth remarking that our model predictions do not 
imply that Lya emitters have low SFRs, but instead have lower 
SFRs than the bulk of the galaxy population at the same UV mag- 
nitude. As shown in the previous section, Lya emitters at high red- 
shifts are predicted to be mainly starbursts. 

Finally, we also study the difference in the size of galaxies. 
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since our outflow geometries depend strongly on the galaxy half- 
mass radius. In the Wind geometry, the hydrogen column density 
scales as ~ and in the Shell geometry the hydrogen 

column density scales as ~ {Ri/2)~^- In order to obtain a high 
Lya escape fraction, Nh has to be low in general, and this could 
be attained by having a large {_Ri/2)- Hence, not surprisingly, our 
outflow geometries predict that Lya emitters typically have larger 
sizes than the bulk of the galaxy population. This result is indepen- 
dent of redshift and UV magnitude, although in detail the difference 
in sizes may v ary with these quantities. A similar result is reported 
bv lBond et alj (|2009j) at z = 3. They find the sizes of Lya emitters 
(in Lya emission) to be always smaller than ~ 3kpc/h, consis- 
tent with the half-mass radius of our typical Lya emitters at this 
redshift. 



7 SUMMARY AND CONCLUSIONS 

In this paper we couple the galaxy properties predicted by the 
iBaugh et all j2005l) version of the semianalytical model GALFORM 
with a Monte Carlo radiative transfer model of the escape of Lya 
photons to study the properties of Lya emitters in a cosmological 
context. 

Motivated by observational evidence that galactic outflows 
shape the asymmetric profiles of the Lya line, we developed two 
different outflow geometries, each defined using the predicted prop- 
erties of galaxies in GALFORM in a slightly different way. Our Shell 
geometry, which consists of an expanding thin spherical shell, has 
a column density Nh proportional to the cold gas mass in the ISM 
of galaxies. Our Wind geometry, on the other hand, consists of a 
spherical expanding wind with number density that decreases with 
increasing radius. The column density in the Wind geometry is re- 
lated to the mass ejection rate from supernovae, which is computed 
by GALFORM. 

We study in detail the Lya line profiles and characterise them 
in terms of their width, asymmetry and offset with respect to the 
line centre, as a function of the outflow column density, expansion 
velocity and metallicity. The Lya properties of the outflow geome- 
tries we study are found to be sensitive to the column density, ex- 
pansion velocity and the geometry of the outflows, as shown in Fig. 
[2I Metallicity is found to have a smaller impact on the line profiles, 
although it has a great impact on the Lya escape fraction, as shown 
in Fig.m 

Both the width and offset from the line centre are found to in- 
crease for outflows with increasing column densities and expansion 
velocities. In both cases, the Lya line profiles in the Shell geome- 
try are more affected by changes of these properties than the Wind 
geometry. 

The Lya escape fraction is found to decrease with metallicity 
and column density, since those two properties are directly propor- 
tional to the optical depth of absorption. Also, higher values of the 
expansion velocity tend to increase the escape fraction. This rather 
complicated interplay between the column density, expansion ve- 
locity, metallicity and geometry justifies our choice for computing 
the Lya escape fraction from a fully-fledged Lya radiative transfer 
model instead of imposing a phenomenological model. 

The drawback is, of course, the difficulty in coupling both 
models in order to obtain a Lya fcsc for each galaxy predicted 
by GALFORM. One critical step is to choose the free parameters 
in our models. In order to set these for each model, we attempt 
to fit the observed cumulative luminosity function of Lya emit- 
ters over the redshift range < z < 7. We find that a single 



choice of parameters applied at all redshifts cannot reproduce the 
observed CLFs. Motivated by the observed difference in sizes be- 
twee n the star forming region s in local and high redshift starbursts 
(e.g. iRuiopakarn et alj l20Ilh . we allow the outflow radii in our 
models to evolve with redshift when the galaxies are starbursts. 
By doing this we can find a suitable combination of parameters 
to match the measured Lya CLF. 

It is worth pointing out that the need to invoke a redshift de- 
pendence in the outflow sizes in starbursts could suggest that other 
important physical processes which determine the escape of Lya 
photons may not be included in this work. Although outflows have 
been proposed in the past as a mechanism to boost the otherwise 
very low escape fraction of Lya photons, other scenarios have been 
proposed. A clumpy ISM could boost the escape of Lya photons 
with respect to that of Lyman continuum photons, since the for- 
mer have a probability of bouncing off a dust cloud, and therefore 
have more chances to escape, whereas in the latter case photons 
travel thr ough dust clouds, and hence have m ore chances of being 
absorbed ( lNeufeld[l99ll ; [Hansen & oi]|2006h . 

Another physical effect not included in our modelling is 
the attenuation of Lya radiation due to the scattering of Lya 
photons when crossing regions of neutral hydrogen in the IGM. 
This effect has been shown to be important at z > 6 when 
the fraction of neutral hydrogen in the IGM is thou ght to 
have been significant (see, e.g. Diikstra. Lidz & Wvithe 2007 
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Laursen. Sommer-Larsen & Razoumovl 2011 ). Hence, we do not 
expect this to significantly change our model predictions at low 
redshifts. 

A more fundamental uncertainty lies in the calculation of the 
intrinsic Lya luminosity of galaxies, which does not depend upon 
the radiative transfer modelling. We can assess whether GALFORM 
computes the correct intrinsic Lya luminosity by studying the in- 
trinsic (unattenuated) Ha luminosity function, since both emission- 
line luminosities are directly related to t he production rate of Ly- 
man continuum photons (see lOrsi et aL I I2OIQ) , and differ only b y 
their case B recombination emission coefficient ( IOsterbrockli98^ . 
Moreover, the Ha emission from galaxies is less sensitive to dust 
than Lya since these photons do not undergo multiple scatterings 
in the ISM, making their path lengths shorter than the typical path 
lengths Lya photons experience, and thus making it easier to esti- 
mate their attenuation by dust. 

The attenuation by dust can be estimated by computing the 
ratio of the intensity of two or more emission lines and comparing 
with the expectation for case B recombination (see, for example, 
iKennicutJ 19831.1 1998h . By compa ring obser ved dust-attenuated Ha 
LFs with GALFORM predictions jOrsi et al . 2010) we have found 
that, at z ~ 0.2, GALFORM roughly overestimates the intrinsic Ha 
luminosities by a factor ~ 3. The scatter in the observed Ha lu- 
minosity functions is large, making it difficult to estimate this ac- 
curately. Nevertheless, the intrinsic Lya luminosities predicted at 
z — 0.2 are a factor ~ 10 or more brighter than what is needed 
to reproduce the observational results, if we fix the free parameters 
of the outflow geometries, meaning that the uncertainty in the in- 
trinsic Lya luminosity is not solely responsible for the discrepancy 
between the observed and predicted LFs. 

Despite this, our simple outflow geometries are found to be in 
agreement with a set of different observations, implying that our 
modelling does reproduce the basic physical conditions determin- 
ing the escape of Lya photons from galaxies. 

A direct prediction of our outflow geometries is the distribu- 
tion of hydrogen column densities of Lya emitters. We find that 
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both models feature column densities with values ranging from 
cm 1 , which is consistent with observational 
estimates shown in I Verhamme et al . (2008). The column density 
distribution is closely related to the predicted Lya escape fractions. 
We find that bright Lya emitters generally have high escape frac- 
tions, and faint Lya emitters have low escape fractions. An im- 
portant consequence of this is that certain galaxy properties which 
correlate with the intrinsic Lya luminosity (such as the instanta- 
neous SFR) do not correlate in a simple way with the observed 
Lya luminosity. In other words, galaxies with the same intrinsic 
Lya luminosity (or, say, the same SFR) could be observed with 
different Lya luminosities, making the interpretation of the prop- 
erties of Lya emitters complicated. 

The predicted escape fractions in the models are remarkably 
consistent with observational measurements (Fig.|9]l, giving further 
support to the scenario of Lya photons escaping through galactic 
outflows. Although the Wind geometry is only partially consistent 
with the observed escape fractions, it is worth pointing out that the 
observational data used to make the comparison in Fig. |9] does not 
constitute a representative sample of the galaxy population at this 
redshift. 

Since our Monte Carlo radiative transfer model records the 
frequency with which photons escape from an outflow, we make 
use of this information to study the predicted Lya line profiles from 
our outflow geometries and compare them with observational mea- 
surements. In some cases the agreement is remarkably good, but 
in others we find significant differences (see Fig. |13t . A detailed 
study of the line profiles could reveal important information about 
the galaxy properties. We plan to undertake such a study of Lya 
line profiles in the context of a galaxy formation model in a future 
paper. 

Finally, our models predict that only a small fraction of galax- 
ies should be selected as Lya emitters. We illustrate in Fig.ll4lthat 
Lya emitters are found in general to have low metallicities (except 
at z = 0.2), low instantaneous SFR and large sizes, compared to 
the overall galaxy population. These constraints arise naturally as 
a consequence of the radiative transfer modelling incorporated in 
GALFORM in this paper. Galaxies need to have low star formation 
rates and large sizes in order to display an outflow with a small col- 
umn density. Since the dust content depends on the optical depth of 
dust and not the metallicity alone, Lya emitters are not necessarily 
low-metallicity galaxies in our models. 

The models presented here for the emission of Lya represent 
an important step towards a detailed understanding of the physi- 
cal properties of these galaxies. With the advent of large observa- 
tional campaigns in the forthcoming years focusing on detecting 
Lya emitters at high redshifts, new data will help us refine and im- 
prove our physical understanding of these galaxies, and thus, enable 
us to improve our knowledge of galaxy formation and evolution, 
particularly in the high redshift Universe. 
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APPENDIX A: DESCRIPTION OF THE MONTE CARLO 
RADIATIVE TRANSFER CODE 

In the context of Lya radiative transfer, photon frequencies, v, are 
usually expressed in terms of Doppler units x, given by Eq. l |20| l. 
The thermal velocity dispersion of the gas, Wth, is given by 

i^th = , (Al) 
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where fcs is the BoUzmann constant, T is the gas temperature, nip 
is the proton mass and is the central frequency of the Lya line, 
va = 2.47 X lO^^Hz. 

When a Lya photon interacts with a hydrogen atom, the scat- 
tering cross section in the rest frame of the atom is given by 



/i 



r/47r2 



mec [u - voY + (r/47r)2 



(A2) 



where /12 = 0.4162 is the Lya oscillator frequency, and T = 
A12 = 6.25 X 10*s~^ is the Einstein coefficient for the Lya tran- 
sition (n = 2 to n = 1). 

The optical depth of a Lya photon with frequency v is deter- 
mined by convolving this cross section with the velocity distribu- 
tion of the gas, 



T,{s) = / / n{V,)a{u,V,) dV,i 

Jo J -00 



(A3) 



where T4 denotes the velocity component along the photon's direc- 
tion. Atoms are assumed to have a Maxwell-Boltzmann velocity 
distribution in the rest frame of the gas. In Doppler units, the opti- 
cal depth can be written as 



r^(s) = aH{x)nHS = 5.868 x 10"-'*r4""'^iVH 



1/2 H{x,a) 



(A4) 



where uh is the hydrogen density, Nh the corresponding hydrogen 
column density, Ta the temperature in units of 10 K and a is the 
Voigt parameter, defined as 



a ^ £Z1^ = 4.7 X 10-V-i/^ 



(A5) 



The Hjerting function H{x,a) jHiertin ell 19381) describes the 
Voigt scattering profile. 



a r^°° e"" dv 



(A6) 



which is often approximated by a central resonant core and power- 
law "damping wings" for frequencies below/above a certain 
boundary frequency Xc, which typically ranges between 2.5 < 
Xc < 4. As a consequence, photons with frequencies close to the 
line centre have a large scattering cross section compared to those 
with frequencies in the wings of the profile. Hence, photons will be 
more likely to escape a medium when they have a frequency away 
from the line centre. 

Scattering events are considered to be coherent (the frequency 
of the photon is the same before and after the scattering event) only 
in the rest frame of the atom, but not in the observer's frame. Thus, 
the thermal motion of the atom, plus any additional bulk motion of 
the gas, will potentially Doppler shift the frequency of the photons, 
giving them the chance to escape from the resonant core. 

In the following, ^i,£,2,^3, ■■■ are different random numbers 
in the range [0, 1]. 

The location of the interaction (with either a dust grain or a 
hydrogen atom) is calculated as follows. The optical depth rint the 
photon will travel is determined by sampling the probability distri- 
bution 



and so 



P{r) 



Tint = -ln(l-Cl) 



(A7) 



(A8) 



This optical depth corresponds to a distance travelled s given by 



r(s) = r^is) + Td{s), 



(A9) 



where Tx (s) and (s) are the optical depths due to hydrogen atoms 
and dust grains respectively. The length of the path travelled is de- 
termined by finding the distance s where r(s) = Tint by setting 



Tint 



(AlO) 



where Ud and ad, the number density of dust grains and cross- 
section for interaction with dust, are described below. 

The new location of the photon corresponds to the point where 
it interacts with either a hydrogen atom or a dust grain. To find out 
which type of interaction the photon experiences, we compute the 
probability Ph (x) of being scattered by a hydrogen atom, given by 



Ph{x) = 



nH(TH{x) 



nHO-H{x) + UdOd 



(All) 



We generate a random number ^2 and compare it to Ph- If ^2 < 
Ph, then the photon interacts with the hydrogen atom, otherwise, 
it interacts with dust. 

When interacting with a dust grain, a Lya photon can be ei- 
ther absorbed or scattered. This depends on the albedo of dust par- 
ticles. At the wavelength of Lya, the albedo is yl ~ 0.4, depend- 
ing on the extinction curve used. If the Lya photon is absorbed, 
then it is lost forever. If not, then it will be scattered. The new di- 
rection will depend on a probability distribution for the elevation 
angle 9, whereas for the azimuthal angle (f> the scattering will be 
uni formly distributed. The scatte ring angle 9 can be obtained from 
the lHenvev & Greensteinl ( Il94lh phase function 



9 



2 (1 + 52 - 2pM)3/2 ■ 



(A 12) 



where — cos 9 and g — (/i) is the asymmetry parameter. If g = 
0, Eg. d A 12b reduces to isotropic scattering, g = 1(— 1) implies 
complete forward (backward) scattering. In general g depends on 
the wavelength. For Lya photons, g — 0.73. 

If the photon is interacting with dust, then we generate a ran- 
dom number ^3 to determine whether it is going to be absorbed or 
scattered, comparing this number to A. If the photon is absorbed, 
then it is lost. If it is scattered, then a new direction must be drawn. 

The interaction of photons with hydrogen atoms is more com- 
plicated. Inside an HI region, atoms move in random directions 
with velocities given by the Maxwell-Boltzmann distribution. Each 
of these atoms will see the same photon moving with a different fre- 
quency, due to the Doppler shift caused by their velocities. Since 
the cross section for scattering depends on the frequency of the 
photon, the probability for an atom to interact with a photon will 
depend on a combination of the frequency of the photon and the 
velocity of the atom. 

To compute the direction of the photon after the scattering 
event, a Lorentz transformation to the frame of the atom must be 
made. The direction of the photon in this frame, n^, is given by a 
dipole distribution, with the symmetry axis defined by the incident 
direction n'- 



P{9) = -(l-fcos^ 
8 



(A13) 



where 9 is the polar angle to the direction n^. The azimuthal angle 
of the outgoing photon is uniformly distributed in the range ^ 
(j) < 2tt. With a Lorentz transformation back to the frame of the 
medium we obtain no, the new direction of the photon. Finally, the 
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new frequency of the photon is given by 

Xf = X — rii ■ u -\- no ■ u, 
= a; — Mil + no ■ u, 



(A14) 
(A 15) 



where u = v /vth is the velocity of the hydrogen atom in units of 
the thermal velocity, and uy is the atom's velocity component along 
the photon's direction prior to the scattering event. 

The algorithm described above will follow the scattering of a 
photon until it escapes (or is absorbed), and then the process starts 
again with a new photon traveling on a different path, and so on un- 
til we are satisfied with the number of photons generated. In prac- 
tice, for the runs shown in this paper the number of photons gener- 
ated varies between a few thousand up to several hundred thousand, 
depending on the accuracy of the result we wish to achieve. 

For the typical HI regions studied here, the number of scat- 
terings that photons undergo before escaping could be as high as 
several tens or hundreds of millions. If we want to model several 
thousand photons, then the total number of calculations grows enor- 
mously and the task becomes computationally infeasible. However, 
most of the scattering events will occur when the photon is at the 
line centre, or very close to it, where the cross section for scattering 
peaks. Whenever the photon falls near the centre it will experience 
so many scatterings that the actual distance travelled between each 
scattering event will be negligible, since in this case it will most 
likely be scattered by an atom with a velocity close to zero. Hence, 
the frequency after such scattering will remain in the resonant core. 
This motivates the possibility of accelerating the code performance 
by skipping suc h inconsequential scattering events. 

Following iDiikstra, Haiman & Spaansl ( [2006) . a critical fre- 
quency, a;crit, defines a transition from the resonant core to the 
wing. Whenever a photon is in the core (with < a;crit) we can 
push it to the wings by allowing the photon to be scattered only by 
a rapidly moving atom. We do this by modifying the distribution 
of perpendicular velocities by a truncated Gaussian, i.e. a distribu- 
tion which is a Gaussian for |m| > Xcrit but zero otherwise. The 
modified perpendicular velocities are then drawn from 



U±2 



ln(e4)sin(27rC2). 



(A16) 
(A17) 



When doing this, we allow the photon to redshift or blueshift away 
from the line centre, thus reducing the cross section for scattering 
and increasing the path length. For the configurations studied here, 
we found that a value of Xcrit = 3 provides a good balance between 
accuracy and efficiency of the code, reducing the execution time by 
a factor 100 or more with respect to the non-accelerated case. 



Al Validation of tlie radiative transfer code 

The flexibility of our Monte Carlo radiative transfer code allows 
us to reproduce configurations for which analytical solutions are 
available. Hence, these are ideal to test the performance and accu- 
racy of the code. In the following we describe the tests we have 
performed on our code, where each comparison with an analytical 
solution tests a different aspect of the code. 

Fig. lAll shows the resulting redistribution function for 3 differ- 
ent initial frequencies using ~ IC photons. There is a remarkably 
good agreem ent between the Monte Carlo code and the analytical 
expression of lHummej j 1 9621) for coherent scattering with a dipolar 
angular distribution, including radiation damping. 

The emergent spectrum from an optically thick, homogeneous 
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Figure Al. The redistribution function of Lya photons scattered by hydro- 
gen atoms for different initial frequencies. The histograms show the result- 
ing frequency distribution from the Monte Carlo code, whereas the dashed 
lines show a numerical integration of the analytical solution of iHummeJ 
0962,). 
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Figure A2. Lyo spectrum emerging from a homogeneous static slab at 
T = 10[K], for optical depths at the line centre of tq = 10^, 10'' and lO'', 
as shown in the plot. The profiles are symmetric around a: = 0. The more 
optically thick the medium, the farther from the line centre the resulting 
peaks of each profile ai'e found. The solid lines show the analytical solution 
by Harrington 1 1973), and the orange and blue histograms show the results 
from the Monte Carlo code for a choice of Xcrit = and 3, respectively. 



static slab with photon s gene rated at the line centre was first cal- 
culated byjHMringtoa l ll973h . and the result was generalised by 
iNeufeld mgoir allowing the generated photons to have any fre- 
quency. 

Fig. IA2I shows the emergent spectrum from a simulated ho- 
mogeneous slab. The temperature of the medium was chosen to be 
T = lOK, since in this regime the analytical expression is accurate 
for optical depths down to tq ^ 10*, which is faster to compute 
with the code. 
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Figure A3. Mean number of scatterings as a function of the optical deptli in 
the line centre of the medium. The circles show the results from the Monte 
Carlo code for config urations wi t h diffe rent ro . The dashed line shows the 
analytical solution o f lHarringtoiJ h973h . 



The typical Lya line profile is double peaked, and is sym- 
metrical with respect to the line centre. The centres of the peaks 
are displaced away from a; = by a value determined by to. The 
higher the optical depth, the farther away from the line centre and 
th e wider the profile w ill be. Fig. IA2 [ compares the analytic solution 
of iHarringtoi] i ll 973b with the ouput from the basic code (orange 
histogram), and the accelerated version (blue histogram). Overall, 
it is clear that the non-accelerated version of the code reproduces 
the analytical formula over the range of optical depths shown here. 
When a;crit = 3 (the blue histograin in Fig. lA2b . the output is vir- 
tually indistinguishable from the non-accelerated version, but the 
running time has been decreased by a factor ~ 200. Therefore, 
Fig. lA2l confirms that the choice of Xcrit = 3 does not compromise 
the accuracy of the res ults. 

iHarringtor] ( Il973h also computed the mean number of scatter- 
ings expected before a Lya photon escapes from an optically thick 
medium for the homogeneous slab. He found 



(Af.cat) = 1.612ro. 



(A 18) 



Fig. I A3 1 shows a comparison between the mean number of scat- 
terings computed using our code with the analytical prediction of 
lHarringtonlll973h . The agreement is remar kably good. 

Fol l owing closely the methodology of Harringtonl ( Il973h and 
iNeufeld ( Il99a) . iDiikstra. Haiman & Spaa irs' feOOe^ computed the 
emergent spectrum from a static sphere. Fig. lA4l shows a compari- 
son between the analytic prediction and the output from the code at 
different optical depths. Again, there is a very good agreement be- 
tween the two. The optical depths shown in Fig |A4l were chosen to 
be different from those in Fig. lA2l to show that the code is follow- 
ing closely the expected emergent spectrum for a range of optical 
dept hs spanning sev eral orders of magnitude. 

iNeufeldId 19901) computed an analytical expression for the es- 
cape fraction of photons emitted from an homogeneous, dusty slab. 

Fig. |A5| shows a comparison between the escape fraction ob- 
tained from a series of simulations, keeping to fixed and varying 
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Figure A4. Lyo spectrum emerging from a homogeneous static sphere at 
T = lOK, for optical depths at the line centre of tq = 10^, 10^ and lO'^. 
The profiles are symmetric around a; = 0. The thicker the medium, the 
farther from the line centre the resulting peak s of each profile are found. 
The so lid lines show the analytical solution of lDiikstra. Haiman & SpaansI 
and the histograms show the results from our Monte Carlo code. 
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Figure AS. The escape fraction of Lya photons from an homogeneous 
dusty slab. The optical depth of hydrogen scatterings at the line centre to 
is held constant at to = 10'^, and different values of the optical depth of 
absorption t^ are chosen. Circles show the output f rom the code, a nd the 
solid orange curve shows the analytical prediction of lNeufelj )l99(ll) . 



Ta, with the analytical solution of iNeufeldl ( Il99 0*). We find a re- 
markable agreement between the analytical solution and our code. 
The escape fraction, as expected, decreases rapidly for increasing 
Ta, which, for a fixed to, translates into having a higher concentra- 
tion of dust in the slab. 

To validate the effect of bulk motions in the gas, we model 
the case of an expanding homogeneous sphere, with a velocity at a 
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Figure A6. The emergent Lya spectrum from a linearly expanding sphere 
with velocity zero at the centre and velocity at the edge Vmax = 
0,20,200 and 2000km/s shown in orange, blue, red and green re- 
spectively. The optical depth at t he line centre is kept fixed at tq = 
-|^q7.06 ■pjjg analytical solution of iDiikstra. Haiman & Spaand {2006) for 
the static case is shown in black. The coloured histograms show the out- 
put from the code. The gray solid curve s show the results obtained with the 
iLaursen. Razoumov & Sommer-LarsenI j2009l) code MoCaLaTA (their Fig. 
8). 



distance r from the centre given by 



VbM = Hr, 

ymax 



H 



R 



(A19) 
(A20) 



where Vmax is the velocity of the sphere at its edge, and R is the 
radius of the sphere. 

There is no analytical s olution for this c onfigur ation (ex- 
cept when T = 0, see iLoeb & RvbickH ( Il999h ). so we 
decided to compare our results to those found by a simi- 
lar Monte Carlo code. We perform this comparison with me- 
dia at T = 10'* [is']. Fig. |A6] shows a comparison between 
our code an d the results obtained with the MoCaLaTA Monte 
Carlo code ( iLaursen. Razoumov & Sommer-LarsenI |2009|) . Both 
codes agree very well. Moreover, the figure helps us understand 
the effect of bulk motions of the gas on the emergent spec- 
trum. First, when Vmax = w e recover the static solution, 
jPiikstra. Haiman & Spaanil2006h . When Vmax ~ 20km/s, the 
velocity of the medium causes photons to have a higher probability 
of being scattered by atoms with velocities dominated by the veloc- 
ity of the medium. These atoms see the photons as being redshifted, 
and hence the peak of the spectrum is shifted slightly towards the 
red part of the spectrum, although still a fraction of photons appear 
to escape blueshifted. When Vmax = 200km/s, the blue peak is 
completely erased, and the peak is shifted even further to the red 
side. For very high velocities, such as Vmax ~ 2000km/s, the 
velocity gradient makes the medium optically thin, and the average 
number of scatterings decreases drastically, and consequentially the 
photons have less chance of being redshifted far into the wings, thus 
shifting the peak back to the centre, but still with no photons in the 
blue side of the spectrum. 



APPENDIX B: THE EFFECT OF THE UV BACKGROUND 

An additional feature of the Wind geometry is the option of com- 
puting the resulting ionization of the medium by photons in the 
intergalactic UV background. The Wind geometry, as described in 
detail in section [2!3.2l is assumed to be completely neutral, but pho- 
toionization from the UV background could have the effect of mod- 
ifying the density profile of the neutral gas. 

It is generally believed that the extragalactic UV background 
is dominated by radiation from quasars and mas s ive young stars 
in sta r forming galaxies l lHaardt & Madaj 1 19961 1200 ll : iMeiksinI 
I2OO9I) . The mean intensity of the UV background at the observed 
frequency vo and redshift 20 is defined as 



Jo(i^o, zo) 



1 

47r 



dz (1 + 2)'^ 

(Bl) 

where z is the redshift of emission, v — voil + z)/{l + zo), 
dl/dz is the line element in a Friedmann cosmology, e is the 
proper space-averaged volume emissivity and Tcff is an effective 
optical depth due to absorption by the IGM. There is no explicit 
solution of equation l IBlb since it must be computed it eratively 
by so lving the cosmological radiative transfer equation jPeeblej 
Il993h. For our analysis w e use the values of Jo (1^,2) tabulated 
bv lHaardt & Madaij JioOlh . Notice, ho wever, that more recent cal- 
culations of the UV background flux (Bolton & Haehnelt 2007; 
Meiks in 2009) show that the p hotoionizing background predicted 
by the lHaardt & Madaul ( l200l[) model may be an underestimate at 
z > 5. 

The fraction of ionized hydrogen x = nHii/nn varies ac- 
cording to a balance between radiative and collisional ionizations 
and recombinations inside the cloud: 



aAUeX = {Th + /3Hne)(l - x), 



(B2) 



where ua = 4.18 x 10 [cm s ] is the c ase A recombination 
coefficient at T = IQ'^K ( IOsterbrock|[l98? ), the photoionization 
rate Th {z) from the UV background is given by 



^h{z) = 



'i-nJo{v, z) 
hv 



<Jv{H)dv, 



(B3) 



and the collisio nal ionizati on rate at T = 10* A", is /?h — 6.22 x 
10"^'^ [cm''s~^]l iCerill 19921) . 

As the UV flux penetrates the outflow, it will be attenuated by 
the outer layers of neutral hydrogen. The UV flux reaching an inner 
layer of the HI region is attenuated by this self-shielding process 
according to 



(B4) 



where Joiy) is the original, un-shielded UV flux, and the optical 
depth t{v) when UV photons travel a distance d inside the HI re- 
gion (coming from outside) is given by 



t{u) = g^{H) 



nH{r)dr. 



(B5) 



flout 



The photoionization rate is computed from the outer radius 
inwards. For each shell inside the outflow, Fh is computed tak- 
ing into account the attenuation given by equations l lB4t and \B5l . 
making the photoionization rate smaller as photons penetrate inside 
the HI region. 

We have found that the result of this calculation modifies 
strongly the outer layers of the wind, but since the number density 
of atoms increase rapidly when going inwards, the self-shielding 
effect effectively suppresses the UV radiation for the inner layers 
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of the wind. For winds with low neutral hydrogen column densities 
the UV background is found to penetrate deeper into the wind, but 
in general the innermost region is left unchanged. 

In principle, we could compute the effect of the UV back- 
ground on the Shell geometry as well. However, we do not per- 
form this calculation since, in this case, the number density inside 
the outflow layers depend strongly on the physical thickness of the 
Shell, which in turn depends on the parameter fth- As discussed in 
the previous section, this parameter is considered to have an arbi- 
trary value provided that fth ^ 0.9. If we include the UV back- 
ground in the Shell geometry, the Lya properties would depend on 
the value of fth assumed, which is an unnecessary complication to 
the model. 



APPENDIX C: A GRID OF CONFIGURATIONS TO 
COMPUTE THE ESCAPE FRACTION 

GALFORM typically generates samples numbering many thousands 
of galaxies brighter than a given flux limit at a number of redshifts. 
The task of running the radiative transfer code for each galaxy is 
infeasible considering the time it takes the Monte Carlo code to 
reach completion, which varies from a few seconds up to several 
hours for some extreme configurations. Therefore, this motivates 
the need to develop an alternative method to assign a Lya escape 
fraction for each galaxy predicted by GALFORM. Instead of running 
the radiative transfer code to each galaxy, we construct a grid of 
configurations spanning the whole range of galaxy properties, as 
predicted by GALFORM. 

The first step to construct the grid is to choose which param- 
eters will be used. In principle, each outflow geometry (Wind or 
Shell) requires 4 input parameters from GALFORM: three of these, 
(K:irc), {R1/2) and (^gas) are used by both geometries. In addi- 
tion, Mej is required in the Wind geometry, and (A/gas) in the Shell 
geometry. 

However, a grid of models using four parameters becomes 
rapidly inefficient when trying to refine the grid. A grid with an 
appropriate binning of each parameter can have as many elements 
as the number of galaxies for which the grid was constructed, and 
hence also becomes prohibitively expensive. 

Therefore, we look for degeneracies in the escape fraction 
when using combinations of the input parameters from GALFORM, 
in order to reduce the number of parameters for the construction 
of the grid. The idea is to find a combination of parameters which, 
when kept fixed while varying its individual components, gives the 
same escape fraction. 

The natur al choice for thi s is to use the column density Nh as 
one parameter. iNeufelj (1 1 99Cf) found that the escape fraction from 
a homogeneous, dusty slab is a function of the optical depth at the 
line centre tq and the optical depth of absorption Ta- Both quanti- 
ties are, in turn, a function of the column density Nh- In the Shell 



geometry, Nh oc Mgas/{i?i/2) , whereas in the Wind geometry 
Nh oc Mej / {{Ri/2)Ve^p) . Although promising, we find that we 
do not recover a constant escape fraction in the Wind geometry 
when the column density is kept fixed while varying its individual 
terms. The reason is that the expansion velocity plays a more com- 
plicated role when computing the escape fraction, with the escape 
fraction increasing rapidly with increasing velocity regardless of 
the other parameters of the medium. 

Therefore, we construct three-dimensional grids for each out- 
flow geometry. We define Cwind = Mcj/Rinn, and Cshdi = 
Afshen/^Jshcii- Then, in the Wind geometry the parameters are 
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Figure CI. Comparison of the escape fraction obtained using a direct cal- 
culation (/osc.rt) or interpolating from a grid (/esc, grid)- Points represent 
a subsample of ~ 1000 galaxies selected from GALFORM at 2 = 0.2. 



Cwind, 14xp and Zgas, whereas in the Shell geometry we use 
Cahoii, Kixp and Zgas. We choose to cover each parameter with a 
bin size appropriate to recover the expected escape fraction with a 
reasonable accuracy when interpolating in the grid, but also ensur- 
ing that the number of grid elements to be computed is significantly 
smaller than the total number of galaxies in the sample. 

We find that, when covering each parameter in logarithmic 
bins of 0.1 we get escape fractions that are accurate enough, and 
the number of elements of the grid we need to compute is usually a 
factor ~ 20 smaller than the total number of galaxies in the sample. 

We fix the number of photons to run for each grid point to 
compute the escape fraction, since this will determine the speed 
with which each configuration will be completed. By studying the 
resulting luminosity function of galaxies, we find that running the 
code with a maximum number of photons Np = 1000 gives results 
which have converged over the range of luminosities observed. This 
means that the minimum escape fraction we are able to compute is 
fcsc = 10""^. Although there are configurations where fcsc can be 
lower than this, they are found not to contribute significantly to the 
luminosity functions. 

Fig. |C II shows an example of the performance of the grid we 
use to compute the escape fraction in the Shell geometry using a 
sub-sample of galaxies from GALFORM aX z = 0.2, chosen in a 
way to cover the entire range of intrinsic Lya luminosities. The 
accuracy of the grid gets progressively worse with lower escape 
fractions, since these have intrinsically larger errors due to the con- 
straint on the maximum number of photons used to compute fosc. 
However, as discussed previously, we found that to reproduce ac- 
curately the luminosity functions there is no need to reduce the size 
of the parameter bins or increase the number of photons used. 
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